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Pref  ace 


The  person  knowledgeable  of  the  general  subject 
of  this  thesis  will  perhaps  notice  a  very  marked 
similarity  to  two  other  theses  accomplished  at  AFIT  by 
prior  students  who  also  had  an  interest  in  the 
correlations  between  the  solar  wind,  the  interplanetary 
magnetic  field  (IMF),  and  geosynchronous  energetic 
electron  count  rates  (fluxes).  These  were  the  work  of 
Captains  Warren  Smith  and  Douglas  McCormick  of  previous 
GSO  classes  83  and  84.  In  this  study,  new  data  were 
gathered  and  analyzed  but,  rather  than  using  the 
classical  methods  of  Smith  and  McCormick,  the  more 
refined  method  of  Time  Series  Analysis  as  described  by 
Box  and  Jenkins  (1976)  was  used  in  an  attempt  to  derive 
models  for  the  energetic  electron  fluxes  (Box  and 
Jenkins,  1976).  It  is  hoped  that  at  least  some  of  the 
models  derived  as  a  result  of  this  thesis  may  be  used  to 
more  successfully  predict  the  level  of  energetic 
electron  activity  with  sufficient  lead  time  to  allow 
efficient  control  of  the  many  satellites  we  have 
operating  in  this  energetic  particle-filled  region  of 
space.  Moreover,  this  thesis  has  been  another  attempt 
to  show  that  some  relationship  exists  between  the  IMF, 
the  soiar  wind,  and  the  energetic  electron  fluxes. 


Perhaps  the  greatest  thing  one  learns  from 
taking  on  a  project  such  as  this  is  not  just  simply  a 
few  facts  about  the  subject  matter.  On  the  contrary, 
what  one  "learns"  for  certain  is  the  virtually 
ingraspable  enormity  of  physical  processes  at  work  in 
the  universe  for  which  God  has  given  us  only  a  few  keys 
to  try  to  understand.  In  short,  I  have  indeed  had  a 
humbling  experience  in  composing  this  study. 

Even  more  humbling  is  the  realization  of  the 
incredible  amount  of  stored  knowledge  that  those  who 
have  helped  me  in  this  endeavor  command  at  their  mere 
fingertips.  Major  Joe  Litko  is  a  brilliant  man  whose 
near  instantaneous  and  extremely  accurate  recollections 
of  the  Box  and  Jenkins  time  series  details  never  ceased 
to  amaze  me.  Concerning  Major  Jim  Lange,  I  can  only 
sing  similar  praises.  His  knowledge  of  the  physics  of 
the  space  environment  is  astounding.  I  am  indeed 
greatly  indebted  for  the  continuing  patience  and  help  of 
these  two  individuals.  I  am  also  obliged  to  give  my 
thanks  to  Mr.  Jim  Ware,  the  AFIT  School  of  Engineering 
computer  consultant  for  his  much  needed  help  in 
transfering  raw  data  from  tapes  to  files  in  my  personal 
computer  accounts.  This  enabled  me  to  comfortably 
analyze  the  data  in  ways  I  am  used  to. 


Anyone  who  nas 
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had  the  misfortune  of  having  to  deal  with  leading  data 
off  a  tape  knows  what  I’m  talking  about.  Mr.  Ware  made 
the  data  reading  process  a  virtual  ease  for  me  with  his 
great  command  of  the  computer  systems  at  AFIT.  Thanks 
also  goes  to  Dr.  Ray  Klebesadel  of  the  Los  Alamos 
National  Laboratory  and  Ms.  Billie  Dolen  and  Mr.  Ralph 
Post  of  the  National  Space  Science  Data  Center  for  their 
very  courteous  and  concerned  help  in  providing  the  raw 
data  for  this  project. 

On  the  home  front,  where  would  I  be  without  the 


comfort  and  restoration  of  my  loving  wife,  Jennifer? 

Even  in  my  darkest  moments  she  was  able  to  help  me  keep 
my  eyes  focused  on  the  goal  and  gave  me  reassurances 
that  I  was  capable  of  completeing  this  project. 

Finally,  thanks  be  to  the  Lord.  May  I  ever  and 
always  be  His  instrument  and  a  part  of  His  plan  for  good 
works . 


Michael  P.  Halpin 
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Abstract 

This  project  used  a  Box  and  Jenkins  time  series 
analysis  of  energetic  electron  fluxes  measured  at 
geosynchronous  orbit  in  an  effort  to  derive  prediction 
models  for  the  flux  in  each  of  five  energy  channels.  In 
addition,  the  technique  of  transfer  function  modeling 
described  by  Box  and  Jenkins  was  used  in  an  attempt  to 
derive  input-output  relationships  between  the  flux 
channels  (viewed  as  the  output)  and  the  solar  wind  speed 
or  interplanetary  magnetic  field  (IMF)  north-south 
component,  Bz,  (viewed  as  the  input).  The  transfer 
function  modeling  was  done  in  order  to  investigate  the 
theoretical  dynamic  relationship  which  is  believed  to 
exist  between  the  solar  wind,  the  IMF  Bz,  and  the 
energetic  electron  flux  in  the  magnetosphere.  The 
models  derived  from  the  transfer  function  techniques 
employed  were  also  intended  to  be  used  in  the  prediction 
of  f 1 ux  va 1 ues . 

The  results  from  this  study  i ndi cat e  that  the 
energetic  electron  flux  changes  in  the  various  channels 


X 


are-  dependent  on  more  than  simply  the  solar  wind  speed 
or  the  IMF  Bz.  Also,  most  of  the  time  series  models 
developed  here  (for  both  the  individual  energetic 
electron  channels  by  themselves  and  those  developed 
through  transfer  functions)  were  not  suitable  for  use  in 
prediction,  since  the  standard  error  of  the  forecasts 
made  using  these  models  was  unacceptably  high.  However, 
a  few  of  the  models  did  merit  possible  consideration  for 
use  in  prediction  of  fluxes.  These  were  the  individual 
time  series  models  for  the  6.6  -  9.7  MeV  channel.  In 
addition,  the  transfer  function  models  developed  using 
the  solar  wind  as  an  input  and  the  6.6  -  9.7  MeV  channel 
as  an  output  may  be  of  possible  use.  The  channel 
containing  electrons  with  energies  between  9.7  -  16  MeV 


was  also  related  to  the  solar  wind  via  a  transfer 
function  with  a  reasonable  forecast  standard  error. 
Finally7,  most  of  the  transfer  function  models  derived 


with  the  solar  wind  considered  as  the  input  to  a  given 
channel  resulted  in  delay  parameters  of  about  2  days 
between  the  input  change  in  solar  wind  velocity  and  the 
observed  output  change  in  electron  flux  which  supports 
findings  from  prior  studies. 
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A  TIME  SERIES  ANALYSIS  OF 


ENERGETIC  ELECTRON  FLUXES  (1.2- 1C  MeY ) 


AT  GEOSYNCHRONOUS  ORBIT 


I.  Introduction 


Background 


In  the  years  immediately  after  World  War  II  and 
the  capture  of  various  German  secret  plans  for  rockets 


from  their  launch  complex  at  Peenemunde  on  the  Baltic 
coast,  the  imagination  of  the  scientific  community  began 
to  be  taken  with  the  possibilities  for  space  exploration 


inherent  in  such  devices  as  the  V-2  which  could  zoom  to 
altitudes  many  miles  higher  than  had  ever  been 
previously  achieved.  Indeed,  Dr.  James  Van  Allen  of  the 
University  of  Iowa  conducted  a  1958  experiment  with  a 
V-2  spinoff,  the  Jupiter,  which  like  its  parent  was 
developed  by  von  Braun  (Boyd,  1974:28).  Van  Allen  put 
together  a  small  package  for  the  Explorer  I  satellite  ir. 
which  he  installed  Geiger  counters  in  the  Jupiter’s 
payload  and  then  helped  launch  it  into  space  to  continue 
the  high  altitude  cosmic  ray  studies  which  he  was  then 
conducting  (Glasstone,  1965:543).  To  Van  Allen’s  surprise, 
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he  found  that  the  term  "space"  apparently  was  a  misnomer  at 
least  in  the  near-earth  sense.  The  Geiger  counters 
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recorded  very  high  amounts  of  energy  which  caused  a 


great  bewilderment  in  the  scientific  community  as  to  the 
energy’s  source.  The?  Van  Allen  belts,  as  they  are  now 
called,  are  regions  of  very  intense  energy  particles 
which  extend  from  altitudes  about  500  miles  above  the 
earth’s  equator  to  over  40,000  miles.  Thus,  man  began 
to  realize  that  this  area  of  "space"  (the  magnetosphere) 
isn’t  just  empty  but  is  literally  filled  with  these 
tiny,  high  energy  particles. 

But  what  are  these  particles,  and  where  do  they 
come  from?  A  good  portion  of  them  are  energetic 
electrons  and  protons  which  have  become  trapped  in  the 
earth’s  magnetosphere.  It  is  believed  that  many  of 
these  particles  are  produced  by  solar  disturbances, 
carried  outward  from  the  sun  by  the  solar  wind,  and  then 
deposited  into  the  magnetosphere  by  a  little  understood 
mechanism  which  allows  the  Interplanetary  Magnetic  Field 
(IMF)  to  "link"  up  with  the  earth’s  magnetic  field 
(Akasofu,  1983:179).  Thus,  the  particles  of  the  solar 
wind  are  allowed  to  enter  the  magnetosphere  where  they 
are  subsequently  accelerated  and  become  "trapped"  along 
the  earth’s  magnetic  field  lines.  The  kinetic  energy  of 
these  particles  is  very  large  which  implies  that  they 
are  extremely  fast  moving.  The  energies  of  the 
particles  which  will  be  dealt  with  in  this  study  are  in 
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the  range  of  1  to  16  MeV.  One  MeV  is  106  electron  volts 
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where  an  electron  volt  is  equal  to  approximately  1.6  x 


10'19  joules  of  kinetic  energy.  Taking  into 


consideration  the  very  small  masses  involved,  it  should 


be  apparent  that  the  velocities  of  these  protons  and 


electrons  are  extremely  high  in  order  to  produce 


energies  of  the  magnitudes  cited  above.  Indeed,  with 


regard  to  the  energetic  electrons  which  will  be  studied 


here,  the  electrons  must  be  considered  in  a  relativistic 


sense,  since  their  velocity  may  approach  and  in  some 


cases  nearly  equal  the  speed  of  light  (3  x  10s 


meters/sec ) . 


The  next  question  concerns  what  the  implications 


are  for  our  operations  in  space  given  all  these 


particles  in  the  near  earth  space  environment.  The 


implications  are  many.  They  can,  for  example,  cause 


difficult  problems  for  our  spacecraft  deployed  in 


geosynchronous  earth  orbits  inside  the  magnetosphere 


High  levels  of  flux  (defined  as  the  measured  number  of 


particles  passing  through  a  unit  area  per  unit  of  time) 


of  these  small  energetic  masses  can  cause  such  maladies 


as:  (1)  spurious  event  sensing  of  surveillance 


satellites,  (2)  uncommanded  tumbling  motions  of  the 


satellite,  (3)  degraded  sensor  and/or  electronics 


capabilities  due  to  particle  energy  absorption  over  a 
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long  period  of  time,  (4)  alterings  of  satellite  surface 
coatings  such  that  the  satellite’s  operating  temperature 
is  increased  above  the  optimum,  and  (5)  radiation 
sickness,  cancer,  or  gene  mutation  for  astronauts 
traveling  for  any  extended  periods  in  such  conditions 
(Lange,  1 982 : I I I -D- 2 1 , 2 2 ;  Spjeldvik  and  Rothwell, 


1983: 113)  . 


Since  it  is  imperative  that  we  have  satellites 


operating  within  the  Van  Allen  belts  in  the 
geosynchronous  region  of  space  (approximately  22,700 
miles  altitude)  for  communications,  surveillance,  and 
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exploration,  and  since  the  exact  mechanism  for  the 
changes  in  flux  of  these  energetic  particles  in  the  Van 
Allen  belts  is  not  well  understood,  statistical  methods 
are  used  in  an  attempt  to  help  explain  particle 
behavior.  Having  some  ability  to  predict  when  changes 
in  flux  levels  might  occur  is  a  very  desirable  goal 
since  such  a  capability  would  allow  (given  sufficient 
lead  time)  more  efficient  control  of  satellites  which 
might  be  exposed  to  potentially  harmful  flux  levels. 

For  example,  in  the  event  of  a  forecast  particle 
substorm,  sensitive  satellite  subsystems  might  be  turned 
off  so  as  to  avoid  the  deleterious  effects  of  increased 
flux  levels  much  as  a  personal  computer  owner  might  turn 
off  his  system  at  the  onset  of  a  thunderstorm  to  avoid 
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the  potential  for  damage  due  to  lightning  strikes. 
Thus,  many  statistical  studies  of  these  energetic 


particle  fluxes  have  focused  on  classical  methods  such 
as  correlation  analysis  or  analysis  of  variance  (AN’OVA) 
in  attempts  to  link  measured  particle  flux  levels 
(recorded  by  a  number  of  satellites  which  the  US 
currently  has  in  geosynchronous  orbit)  with  other 
magnetospheri c  parameters  which  are  also  measurable  by 
our  satellites  or  at  ground  stations.  The  idea  is  that 
if  a  given  magnetospheric  or  solar  wind  parameter  can  be 
used  with  sufficient  lead  time  (say  one  to  two  days)  to 
predict  an  increase  in  particle  flux,  then  the  parameter 
may  be  used  as  a  basis  for  satellite  command  and  control 
dec i s ions . 

One  of  the  key  areas  of  interest  has  been  the 
belief  on  the  part  of  many  geophysicists  that  there  is  a 
relationship  between  changes  in  the  solar  wind  velocity, 
the  north-south  component  of  the  solar  wind  (or 
interplanetary)  magnetic  field  (IMF),  and  the  levels  of 
flux  of  energetic  electrons  in  the  magnetosphere 
(Paulikas  and  Blake,  1978:2,10).  Studies  of  note  which 
have  made  attempts  to  correlate  such  parameters  via  the 
classical  methods  mentioned  above  are  those  of  Smith 
(1983)  and  McCormick  (1984).  Smith,  employing  data  from 
June  1979  to  April  1982,  used  a  correlation  analysis  and 
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ANOYA  of  solar  wind  speed,  IMF,  and  energetic  electron 
flux  level  data  and  concluded  that  there  was  only  a  weak 
correlation  (largest  value  of  R2  =  0.2582  when  two-day 
old  solar  wind  data  used)  between  solar  wind  speed 
readings  and  energetic  electron  flux  levels  (Smith, 
1983:42-50).  Also,  no  significant  correlations  between 
the  IMF  north-south  (Bz)  component  and  e~  flux  levels 
were  noted  (Smith,  1983:54).  McCormick,  with  the  same 
data  set,  used  methods  similar  to  those  of  Smith  and 
came  to  virtually  identical  conclusions,  but  he  also 
found  that  solar  wind  correlations  with  increased  flux 
levels  were  slightly  higher  when  the  IMF  north-south 
component  was  negative  (ie,  more  southward  oriented). 

To  elaborate,  R2  values  were  0.269  and  0.313 
respectively  for  the  correlations  between  squared  two 
day  old  solar  wind  speed  values  and  energetic  electron 
flux-  values  from  the  two  lowest  electron  flux  channels 
of  1.2  -  1.8  MeV  and  3.4  -  4.9  MeV.  These  R2  values 
were  attained  only  when  the  IMF  north-south  (Bz) 
component  was  negative  during  the  same  time  at  which  the 
solar  wind  speed  was  measured  (McCormick,  1984:51-63). 
This  tended  to  support  existing  ideas  that  a  southward 
oriented  Bz  component  of  the  IMF  seems  to  result  in  more 


enhanced  flux  levels  (Potemra,  1983:279). 


Despite  these  sometimes  discouraging  results, 
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the  scientific  community  continues  to  believe  that  some 
sort  of  relationship  exists  between  the  IMF  Bz,  the 
solar  wind,  and  the  energetic  electron  fluxes,  and 
further  study  in  this  area  seems  warranted  (Paulikas  and 
Blake,  1978:1).  In  fact,  a  number  of  other  studies  have 
very  much  established  that  a  statistical  relationship 
does  exist  when  electron  fluxes  of  lower  energy  are 
involved  (Su  and  Konradi,  1979:25). 


Objectives  of  the  Research 

A  comparatively  new  and  unused  method  of 
statistical  analysis  for  data  taken  over  a  long  time 
period  is  that  of  Time  Series  Analysis  such  as  that 
described  by  Box  and  Jenkins  (1976).  With  the  advent  of 
very  powerful  computers  and  the  availability  of 
statistical  programs  such  as  the  BMDP  Statistical 
Software  package,  the  capability  to  do  analyses  like 
that  of  Box  and  Jenkins  now  exists.  Therefore,  the 
objectives  of  this  research  will  be  twofold.  First,  a 
Box  and  Jenkins  time  series  analysis  of  energetic 
electron  flux  data  will  be  performed  in  an  attempt  to 
de fine  model s  which  may  be  used  to  predict  changes  i n 
flux  levels .  As  stated  earlier,  such  models  might  allow 
the  prediction  of  potentially  dangerous  fluxes  with 


sufficient  lead  time  to  allow  better  control  of  affected 
spacecraft.  Second,  since  one  of  the  main  concerns  of  the 
geophysical /as t rophys i ca 1  community  is  in  the  linkage 
between  the  solar  wind  velocity,  the  north-south  (Bz) 
component  of  the  IMF,  and  the  changes  in  flux  levels  of 
energetic  e-,  an  attempt  will  also  be  made  in  this  study 
to  relate  these  entities  via  time  series  transf er 
functions .  The  scope  of  this  study  will  be  limited  to 
data  taken  in  the  time  period  between  April  1982  and  May 
1986,  a  period  somewhat  longer  than  that  in  the  studies  by 
Smith  and  McCormick.  It  should  be  emphasized  that  some  of 
this  data  may  be  discarded  due  to  the  way  in  which  a  time 
series  analysis  is  conducted.  The  finer  details  of  a  time 
series  analysis  as  described  by  Box  and  Jenkins  will  be 
presented  in  the  Detailed  Methodology  section  of  this 
report . 

Overview 

This  report  includes  a  literature  review  of 
articles  on  geosynchronous  particle  flux  and  its 
implications  (Chapter  II),  data  preparation  and  use  of 
computers  (Chapter  III),  a  detailed  look  at  the  Box  and 
Jenkins  time  series  analysis  (Chapter  IV),  a  presentation 
of  results  (Chapter  V),  and  finally,  appropriate 
conclusions  and  recommendations  (Chapter  VI). 


II.  Literature  Review 
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One  of  the  most  obvious  points  concerning  a 
literature  review  in  the  general  field  of  magnetospheric 
phenomena  is  the  abundance  of  articles.  This  is  due  in 
part  to  the  commissioning  in  the  mid-1970s  of  the 
International  Magnetospheric  Survey  (IMS)  "in  which  a 
coordinated  effort  was  made  to  understand  magnetospheric 
processes"  (Russell  and  Southwood,  1982 :vii).  The  IMS 
was  created  by  a  group  of  concerned  geophysicists  who 
felt  that  more  research  was  needed  concerning  the 
magnetosphere.  The  coordination  of  the  many  projects 
resulting  from  the  multinational  commitment  to  the  IMS 
has  been  and  continues  to  be  carried  out  at  the  level  of 
the  participating  scientists  from  each  nation.  The 
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result  of  the  IMS  has  been  a  virtual  (and  much  needed) 
flood  of  research  projects  and  interest  in  the 
magnetosphere  and  its  mechanics. 

One  of  the  best  publications  for  finding 
literature  on  magnetospheric  processes  is  the  space 
physics  portion  (blue  colored  volumes)  of  the  Journal  of 
Geophys i cal  Research .  This  is  a  monthly  publication 
which  is  filled  with  the  s tudi es  of  the  world’s  leading 
space  physicists.  Many  of  the  articles  cited  here  came 
from  this  journal. 


The  way  in  which  one  chooses  to  conduct  a 


literature  review  of  the  magnetosphere  is  dependent  on 
the  specific  project  concerns.  In  this  case,  it  was 
decided  to  look  at  articles  and  books  associated  with 
the  general  theme  of  the  interactions  between  the  IMF, 
the  solar  wind,  and  the  flux  levels  of  energetic 
particles  observed  in  the  magnetosphere.  Also,  were  it 
not  for  the  fact  that  these  energetic  particles  pose 
difficult  problems  for  our  spacecraft,  one  might  easily 
surmise  that  the  intensity  of  research  in  this  area 
would  be  somewhat  lessened.  Consequently,  reviews  of 
articles  concerning  the  difficulties  of  space  operations 
in  the  magnetosphere  along  with  ways  of  predicting 
hazardous  conditions  there  also  seem  appropriate. 

Magnet ospher ic  Interactions 

A  generally  accepted  theory  is  that  the  particle 
energy  flux  in  the  magnetosphere  is  greatly  controlled 
by  an  as  yet  undetailed  mechanism  whereby  the 
interplanetary  magnetic  field  lines  link  with  the 
earth’s  magnetic  field  lines.  According  to  Nishida, 
extensive  examinations  of  both  ground  based  and  space 
based  observations  have  led  us  to  believe  that  the 
energy  supplied  to  the  magnetosphere  "proceeds  mainly  by 
reconnection  between  the  lines  of  force  of  the  IMF  and 
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the  geomagnetic  field"  (Nishida,  1983:185).  When  this 
occurs,  vast  amounts  of  energy  in  the  form  of  particles 
carried  along  by  the  solar  wind  are  allowed  to  enter  the 
magnetosphere  where  subsequently  these  particles  are 
accelerated  to  extreme  speeds  and  become  trapped  along 
the  earth’s  field  lines.  Pauiikas  and  Blake,  in  an 
earlier  study  of  11  years  of  data  on  electron  fluxes  at 
geosynchronous  orbit,  found  that  the  "efficiency  of 
coupling"  between  these  solar  wind  particles  and  the 
magnetosphere  is  apparently  controlled  by  the  "IMF 
direction  as  well  as  the  solar  wind  velocity”  (Pauiikas 
and  Blake,  1978:2). 

Part  of  the  problem  in  developing  a  deeper 
understanding  of  the  processes  controlling  the  rate  of 
energy  inflow  to  the  magnetosphere  has  to  do  with  the 
vastness  of  space  itself.  According  to  Baker,  the  very 
great  distances  involved  make  it  extremely  difficult  to 
"probe  the  system  concurrently  at  enough  different 
points  to  truly  understand  the  complex  relationships 
between  its  different  parts"  (Baker,  1982b:5917).  Su 
and  Konradi ,  in  an  earlier  paper,  were  in  agreement  and 
stated  that  "the  observations  made  by  a  single 
spacecraft  so  far  fail  to  resolve  the  temporal  and 
spatial  variations  of  the  environment"  ( Su  and  Konradi, 
1979:23).  Nevertheless,  by  correlating  data  from 


1982b:5917).  Nishida  again  echoed  this  feeling  in  a 
later  report  (Nishida,  1983:185). 

Still  another  facet  in  understanding  the 
interactions  is  the  part  played  by  the  north-south  (Bz) 
component  of  the  IMF.  Numerous  studies  have  shown  that 
increases  in  energetic  particle  flux  seem  to  correlate 
well  with  a  negative  (or  southward  oriented)  Bz. 
Researchers  involved  in  such  studies  included  McPherron 
who  in  observations  of  substorms  foui.d  that  the  growth 
phase  of  these  storms  occurred  during  southward  oriented 
Bz  components  of  the  measured  IMF  (McPherron  et  al , 
1973:3131).  Studies  conducted  by  Russell  in  the 
following  year  came  to  the  same  conclusions  (Russell  et 
al ,  1974:1108).  Moreover,  a  1977  study  by  Caan  and 
others  again  reached  the  same  conclusions  regarding  the 
onset  of  substorms  following  a  prolonged  (two  hours  or 
more)  southward  turning  of  the  IMF  Bz  (Caan  et  al , 
1977:4837).  More  recently,  statements  by  Akasofu  and 
Baker  on  this  subject  only  serve  to  emphasize  the 
unanimity  of  the  geophysical  community  concerning  the 
appropriateness  of  this  theory.  In  a  1983  paper, 

Akasofu  wrote, 

.  .  .  it  has  been  found  that  the  north-south 

component  of  the  solar  wind  (or  interplanetary) 
magnetic  field  is  one  of  the  most  important 
parameters  which  link  solar  wind  disturbances 
and  magnetospher  ic  disturbances.  .  .  Thus,  it  is 


crucial  to  understand  the  physical  causes  of 
theta  (the  north-south  component  of  the  solar- 
wind  magnetic  field)  and  its  time  variations  in 
linking  solar  wind  disturbances  and 
magnetospheric  disturbances  (Akasofu, 

1983 : 179 , 181  )  . 

Baker  was  more  to  the  point  when  in  a  later  article  that 
same  year  he  stated, 

.  .  .  a  southward  IMF  (B2  <  0)  gives  rise  to 

enhanced  magnetospheric  activity,  while  a 
prolonged  northward  IMF  (B2  >  0)  is  followed  by 
quiet  geomagnetic  conditions  (Baker  and  others, 
1983:6230) . 

There  are  many  more  geophysicists  in  the  available 
articles  who  are  of  the  same  beliefs. 

Another  point  is  also  brought  out  in  the 
literature  and  bears  mentioning  here.  In  a  study  to 
determine  the  effects  of  the  solar  wind  on 
magnetospheric  dynamics,  Paulikas  and  Blake  noted  that 
increases  in  the  solar  wind  velocity  as  well  as  the 
aforementioned  southward  oriented  IMF  correlated  with 
increases  in  the  numbers  of  energetic  electrons  in  the 
magnetosphere  (Paulikas  and  Blake,  1978:2).  In  fact, 
they  were  so  adamant  regarding  the  significance  of  the 


solar  wind  velocity  that  they  went  on  to  state  that  the 
"velocity  of  the  solar  wind  is  the  most  important 
parameter  in  organizing  the  flux  levels  of  energetic- 
electrons  in  the  outer  magnetosphere"  (Paulikas  and 
Blake,  1978:16).  As  part  of  their  study  they  also 


showed  that  increases  in  the  flux  levels  of  energetic 
electrons  (  >  3.9  MeV  )  correlated  well  with 

approximately  2  day  old  solar  wind  speed  data. 
Specifically, 

Starting  about  a  day  or  two  after  the  solar  wind 
stream  first  reaches  the  earth,  the  fluxes  build 
up  in  correlation  with  the  increasing  velocity 
of  the  solar  wind  (Paulikas  and  Blake,  1978:11). 

This  is  in  general  agreement  with  the  findings  of  Smith 

and  McCormick  who,  as  noted  in  Chapter  1,  found  some 

amount  of  correlation  between  two  day  old  solar  wind 

speed  data  and  increased  flux  level s .( Smi th ,  1983  and 

McCormick,  1984).  One  interesting  side  light  to  the 

Paulikas  and  Blake  paper  has  to  do  with  their 

observation  that  the  most  marked  variability  of  higher 

energy  electron  fluxes  is  associated  with  the  27  day 

solar  rotation  period  (Paulikas  and  Blake,  1978:22). 

This  is  a  characteristic  which  will  clearly  present 

itself  when  the  results  of  the  analysis  are  presented 

later  on  in  this  report. 

Thus  one  can  summarize  the  literature  that 
surrounds  the  interaction  aspect  with  three 
observations:  (1)  the  solar  wind  is  the  energy  provider 
(2)  the  IMF  linkage  with  the  geomagnetic  field  is  the 
key  which  apparently  turns  the  energy  flow  on  and  off, 
and  (3)  the  solar  wind  velocity,  perhaps  more  than  any 
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the  major  concerns  is  "the  on-orbit  lifetime  of 


microelectronic  devices  that  are  designed  to  a  specific 
level  of  radiation  ’hardness’  (such  as  104-10*  rad)" 
which  dictates  a  "trade-off  between  orbit  choice  and 
system  lifetime"  (Spjeldvik  and  Rothwell,  1983:114). 
Since  just  a  single  energetic  electron  of  the  type 
analyzed  in  this  study  may  carry  as  much  as  16  MeV ,  it 
is  easy  to  see  why  space  operations  in  the  magnetosphere 
can  not  be  taken  lightly.  In  addition,  the  study  noted 
that  the  upper  limit  for  human  tolerance  to  such 
radiation  is  only  about  500  rads  (lethal)  while  lesser 
amounts  can  cause  serious  biological  damage  (such  as 
gene  itation  or  cancer)  (Spjeldvik  and  Rothwell, 

1983 :  121  )  . 

Therefore,  given  that  man  has  found  it  necessary 
to  operate  spacecraft  in  the  magnetosphere  for  a  host  of 
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reasons,  it  will  behoove  us  all  to  understand  as  much 
about  magnetospheric  processes  as  possible  and  to 
develop  relationships  for  the  prediction  of  such  events 
so  that  sufficient  lead  time  is  available  for  spacecraft 
control.  To  a  great  extent,  the  practical  approach  up 
till  now  has  been  the  observation  of  significant 
magnetospheric  events  and  the  correlation  of  these  with 
measured  geophysical  parameters. 
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Prediction  Methods 

Numerous  studies  have  been  conducted  for 
determining  the  correlations  between  events  observed  in 
the  magnetosphere  and  measured  parameters.  An  example 
of  one  of  these  studies  attempted  to  prove  or  disprove 
previous  correlations  between  such  things  as  observed 
spacecraft  anomalies  and  local  time,  geomagnetic 
activity,  the  day  of  the  week,  and  the  season  (Grajek 
and  McPherson,  1979).  The  study  concluded  that  the 
anomalies  are,  with  99.997%  confidence,  dependent  upon 
geomagnetic  activity  as  measured  by  the  index  Dsi  where 
Dst  is  a  measure  of  the  equatorial  disturbance  produced 
by  magnetic  storms  (Grajek  and  McPherson,  1979:774).  Su 
and  Konradi  also  derived  a  model  for  particle  flux 
intensities  at  geosynchronous  altitude  using  a  thirl 
order  polynomial  curve  of  best  fit  for  their  available 
data  concerning  particles  of  the  energy  range  70  - 
41,000  eV  (Su  and  Konradi,  1979:27).  During  their  study 
they  also  showed  that  geosynchronous  particle  flux 
intensities  correlated  well  enough  with  the  auroral 
electrojet  (AE)t  index  to  indicate  a  definite  causal 
relationship,  while  the  flux  intensities  did  not 
correlate  well  with  the  geomagnetic  planetary  (Kp)  index 
(Su  and  Konradi,  1979:26,27). 


Others  have  taken  the  approach  of  trying  to 
directly  include  characteristics  of  the  solar  wind  in 
their  representations  of  magnetospher i c  processes.  In  a 
1983  paper,  Akasofu  derived  an  equation  for  total  energy 
output  rate  of  the  magnetosphere  which  contained,  among 
other  factors,  the  solar  wind  velocity  and  the  square  of 
the  overall  solar  wind  magnetic  field  magnitude 
(Akasofu,  1983:176).  McCormick  pointed  out  that  this 
approach  considered  the  magnetosphere  to  be  a  "driven" 
system  rather  than  an  "unloading  system  where  substorms 
occur  from  energy  accumulated  in  the  magnetosphere  and 
(are)  released  by  some  instability"  (McCormick, 

1984:16).  Crooker  and  others  after  looking  at  six-month 
and  yearly  solar  wind  speed  averages  found  a  high 
correlation  with  geomagnetic  activity  and  subsequently 
suggested  that  the  product  of  the  southward  component  of 
the  IMF  and  the  square  (or  higher  power)  of  the  solar 
wind  velocity  seemed  to  correlate  well  with  geomagnetic 
activity  (Crooker  et  al,  1977:1933-1936). 

Alas,  there  are  many  theories  concerning  ways  to 
predict  particle  flux  in  the  magnetosphere,  but  as  yet 
there  is  no  universally  acceptable  way  to  predict,  much 
less  describe,  magnetospher ic  processes.  Perhaps 
Paulikas  and  Blake  summed  up  the  hopes  of  many  space 
researchers  when  they  stated: 


.  .  .  it  seems  clear  that  the  present  results 

already  offer  some  hope  of  both  short-term  and 
long-term  prediction  of  the  energetic  electron 
radiation  in  the  outer  magnetosphere  if  a 
sufficiently  accurate  prediction  of  the 
parameters  of  the  solar  wind  are  available 
(Paulikas  and  Blake,  1978:33). 

With  the  hope  of  deriving  prediction  models  for 

energetic  electron  flux  in  the  magnetosphere  based  on 

the  Box  and  Jenkins  time  series  analysis  methods  along 

with  the  same  hope  for  deriving  a  relationship  between 

the  flux  and  the  IMF/solar  wind  speed  via  transfer 

functions,  let  us  go  on  to  consider  how  the  data  were 


prepared . 


Data ,  Software ,  and  Computers 


Data 

The  data  for  this  project  were  supplied  by  two 
sources.  The  Los  Alamos  National  Laboratory  and  the 
satellite  1979-053  were  the  source  of  one  set  of  the 
data.  The  National  Space  Science  Data  Center  through 
readings  from  the  IMP  8  and  ISEE  3  satellites  was  the 
other . 

The  Los  Alamos  data  was  originally  sent  on  a 
magnetic  tape  and  contained,  among  ten  variables,  daily 
average  values  of  f ive  energetic  electron  flux  level 
channels  covering  the  time  period  from  April  1982  to  May 
1986.  The  details  of  the  orientation  of  satellite 
1979-053  and  its  sensors  which  record  the  flux  values 
are  available  from  Baker  et  al .  The  five  channels 
mentioned  contain  measurements  for  flux  levels  (count 
rates)  in  the  ranges  of  1.2  -  1.8  Mev,  3.4  -  4.9  MeV, 

4.9  -  6.6  MeV,  6,6  -  9.7  MeV,  and  9.7  -  16  MeV.  Two 
different  detector  packages  aboard  the  satellite 
collected  this  data.  The  lowest  energy  channel  (1.2  - 
1.8  MeV)  was  collected  by  a  solid  state  detector. 
Henceforth,  this  channel’s  data  shall  be  referred  to  as 
the  SEESSD  data  channel  (where  "SEE"  is  an  acronym  for 
spectrometer  for  extended  electron  measurements  and 


"SSD"  stands  for  solid  state  detector)  (Baker  et  al  , 
1982a:83).  The  last  four  channels  were  collected  by  the 
other  onboard  detector  and  will  be  referred  to 
henceforth  as  the  SEEI  ,  SEEII,  SEEIII,  and  SEE1V  data 
channels  respectively.  One  point  which  should  be  made 
is  that  the  flux  values  measured  are  simply  the  count 
rates  which  are  not  the  same  as  the  standard  units  for 
flux  levels  (particles  per  unit  area  per  unit  time). 
However,  the  correlation  between  count  rates  and  the 
standard  units  is  direct,  so  the  measured  count  rates 
may  be  used  in  an  analysis  just  as  the  data  converted  to 
standard  units  could  be  used. 

Regarding  the  second  set  of  data  supplied  by  the 
National  Space  Science  Data  Center  (NSSDC),  once  again, 
the  data  were  sent  on  a  magnetic  tape  referred  to 
commonly  by  the  NSSDC  as  the  "Omni  Tape".  This  tape,  in 
addition  to  containing  the  pertinent  values  for  solar 
wind  speed  and  the  IMF  components,  contained  readings  of 
no  less  than  37  different  parameters.  The  big 
difference  between  this  data  and  that  supplied  by  Los 
Alamos,  however,  was  that  the  Omni  Tape  contained  hourly 
readings  of  all  these  different  parameters  over  a  period 
stretching  from  April  1981  to  April  1985.  This  meant 
that  the  data  in  its  raw  form  was  contained  in  a  five 
megabyte  file  with  much  of  the  information  irrelevant  to 
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this  study.  Consequently,  a  computer  program  had  to  be 
written  to  allow  compression  of  the  original  Omni  Tape 
supplied  into  a  more  manageable  file  containing  daily 


averages.  This  was  done  in  order  to  synchronize  daily 
average  values  for  the  solar  wind  speed  and  IMF  with  the 
daily  average  readings  contained  on  the  Los  Alamos  tape 
for  the  energetic  electron  fluxes.  Though  the  NSSDC 
obtains  their  Omni  Tape  data  from  a  total  of  17 
different  satellites,  the  majority  of  the  readings  for 
the  solar  wind  speed  and  the  IMF  component  values  come 
from  the  IMP  8  (for  International  Magnetospheric  Probe) 
and  ISEE  3  (for  International  Sun-Earth  Explorer) 
spacecraft.  The  ISEE  project  is  a  joint  NASA-European 
Space  Agency  effort  to  study  the  outer  magnetosphere 
(von  Rosenvinge,  1982:1).  Once  again,  details 
concerning  spacecraft  orbital  parameters  and 
instrumentation  are  available  from  McCormick  (McCormick, 
1984:7)  and  King  (King,  1982:10-20). 


Overview 


As  a  first  action,  both  tapes  were  read  into 
personal  files  on  the  CDC  Cyber  6000  and  SSC  Unix  VAX 
11/780  computers  at  AFIT.  Next,  due  to  the  intractable 
nature  of  the  data  contained  on  the  Omni  Tape,  its  data 
was  compressed  into  daily  averages  and  then  saved  in  a 
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dif ferent  file. 


The  Los  Alamos  data  for  the  energetic 


electron  fluxes  (count  rates)  was  already  in  this  daily 
form.  A  separate  file  containing  time  synchronized 
solar  wind  speed,  IMF  B2,  and  energetic  electron  flux 
values  was  also  created  to  allow  transfer  function 
modeling.  Once  this  was  accomplished,  analyses  of  the 
data  began.  The  BMDP  Statistical  Software  package  used 
in  this  study  was  available  on  both  systems.  Also, 
files  were  transferable  between  the  two  computers.  The 
availability  of  both  systems  for  this  project  helped 
reduce  time  delays  when  one  system  was  busy.  This  was 
helpful  in  that  many  analyses  were  necessary  with  a 
great  deal  of  iterative  interactions  on  the  part  of  the 
author . 

Detai led  Preparations 

As  explained  earlier,  the  data  for  this  study 
were  obtained  from  LANL  and  NSSDC  in  the  form  of  tapes. 
The  LANL  tape  contained  the  data  for  daily  average  count 
rates  in  five  different  energetic  electron  energy  ranges 
(or  channels).  This  tape  was  obtained  first  and  was 
delivered  to  the  AFIT  school  of  engineering  computer 
consultant  who  was  able  to  successfully  load  the  tape  on 
the  school’s  tape  drive  unit  in  the  computer  terminal 
room.  Using  a  program  to  transfer  data  from  a  tape  file 
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to  a  disk  file  in  the  SSC's  memory,  the  consultant  then 
transfered  the  data  to  a  personal  file  on  the  SSC  which 
was  subsequently  copied  by  the  author  to  enable  later 
analysis  via  BMDP .  The  data  supplied  by  LANL  are 
explained  in  more  detail  in  an  article  by  Baker  (Baker, 
1982a :82-90 ) . 

The  data  supplied  by  the  NSSDC  (Omni  Tape)  also 
came  in  magnetic  tape  form.  The  problem  with  this  file 
was  its  great  size  (approximately  5.2  megabytes).  The 
process  followed  for  "getting"  the  data  from  this  tape 
was  somewhat  different  from  that  followed  for  the  LANL 
tape.  The  tape  was  loaded  on  the  school’s  tape  drive 
unit  in  the  computer  room.  However,  since  its  size  was 
too  great  to  be  accommodated  in  a  personal  file,  the 
consultant  loaded  it  into  a  general  purpose  directory  of 
the  SSC’s  memory.  This  directory  is  purged  every  48 
hours,  but  its  availability  allowed  the  author  enough 
time  to  write  a  small  FORTRAN  program  which  could  access 
the  data  from  the  Omni  Tape  (as  written  in  the  general 
purpose  directory)  and  then  compress  it  into  daily  (24 


rV 


hour)  averages.  Once  this  was  done,  the  size  of  the  new 
file  was  reduced  by  1/24,  and  this  made  it  manageable 
and  small  enough  to  be  stored  in  the  author’s  personal 
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SSC  account.  All  daily  average  parameters  (from  the 
original  37)  not  seen  to  be  of  any  use  to  this  study 


were  discarded  in  order  to  further  reduce  the  size  of 


the  file.  For  completeness  program  readom ,  a  listing  of 
the  FORTRAN’  program  written  to  compress  the  Omni  Tape, 
is  included  in  Appendix  A  of  this  report.  Adequate 
comments  are  available  within  the  program  listing  to 
enable  the  reader  to  understand  what  was  done  to 
actually  read  and  compress  the  NSSDC  data.  Details  of 
the  data  contained  on  the  Omni  Tape  are  contained  in  a 
separate  set  of  articles  (von  Rosenvinge,  1982:1-9  and 
King,  1982:11-20). 

A  set  of  particulars  concerning  each  tape’s 
generation  also  accompanied  them  in  the  mail.  This  was 
both  helpful  and  necessary  in  each  case,  since  without 
knowledge  of  the  format  of  the  information  on  each  tape, 


the  writing  of  a  program  to  read  the  data  would  have 
been  impossible. 

In  order  to  perform  a  time  series  analysis  using 
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transfer  functions ,  it  was  necessary  to  write  another 


FORTRAN  program  which  could  read  values  from  the  NSSDC 


file  (solar  wind/IMF  values)  and  the  LANL  file  (electron 


flux  values)  and  then  store  these  values  in  a  single 


time  synchronized  data  file.  This  was  required  so  that 


computation  of  the  cross  correlation  f unct i ons  required 


in  transfer  function  modeling  could  take  place.  Only 


the  values  from  overlapping  days  of  the  two  files  were 


ism 


read  and  stored,  since  the  cross  correia* ion  between  two 
time  series  only  makes  sense  if  the  series  are 
synchronized  in  time.  Thus,  progam  combo  1  read  and 
stored  706  combined  cases  of  data  from  May  8,  1983  to 
April  12,  1985.  This  encompassed  approximately  the  last 
two  years  of  the  NSSDC  data  on  solar  wind  speed  and  IMF 
and  the  middle  two  years  of  the  LANL  data  on  the 
energetic  electron  flux.  This  was  more  than  sufficient 
to  perform  a  suitable  transfer  function  analysis.  A 
copy  of  program  combol  is  also  included  in  Appendix  B. 

Missing  Values 

Some  important  points  concerning  the  data  should 
be  mentioned  which  have  a  direct  bearing  on  this  study. 
On  both  tapes  some  data  were  naturally  missing.  This  is 
of  course  due  to  the  inherent  problems  with  trying  to 
obtain  satellite  readings  via  sensors  which  can 
malfunction.  Missing  data  can  have  serious  implications 
for  a  time  series  analysis  which  may  only  be 
accomplished  on  a  set  of  equally  spaced  (in  time) 
readings.  To  remedy  this  problem,  two  different  methods 
were  used.  Regarding  the  energetic  electron  data  (the 
tape  from  LANL),  missing  values  were  estimated  by  BMDP 
program  AM,  "Description  and  Estimation  of  Missing  Data" 
using  the  SINGLE  method  whereby  'a  value  for  a  variable 


is  estimated  by  regressing  that  variable  on  the  variable 


with  which  it  is  most  highly  correlated"  (Dixon  et  al , 
1985:217-234).  This  seemed  like  a  logical  method  since 
the  data  contained  five  different  electron  flux  channels 


with  at  least  some  amount  of  presumable  correlation 
between  the  channels.  The  use  of  program  AM  was 
possible  because  there  was  a  total  of  only  32  days  over 


the  four  year  period  covered  by  the  data  where  missing 
values  occurred.  Also,  the  smaller  size  of  the  LANL 
file  made  it  easy  to  apply  program  AM.  The  majority  of 


these  missing  values  were  in  August  of  1982  (17  were 
missing)  which  was  subsequent  cause  for  disregarding  the 
data  from  1982  for  this  study.  The  estimated  values 


provided  by  program  AM  were  then  substituted  into  the 
data  set  to  allow  analysis. 

The  NSSDC  data,  however,  presented  a  slightly 


different  problem.  The  massive  size  of  the  Omni  Tape 
made  its  data  intractable  as  far  as  using  BMDP  program 
AM.  Therefore,  the  author  decided  that  the  best  way  to 
handle  missing  data  on  the  Omni  Tape  was  to  allow 
program  readom  to  store  the  last  computed  daily  average 


as  the  new  daily  average  if  an  entire  day’s  set  of 
hourly  readings  was  missing.  Also,  if  only  part  of  one 
day’s  hourly  readings  was  available  for  averaging  ( say 
only  10  hourly  readings  instead  of  the  usual  24  for  one 


day),  then  those  readings  which  were  available  were 


averaged  and  taken  to  be  the  daily  average  value. 

Again,  only  in  the  cas<.  where  an  entire  day's  set  of 
hourly  readings  was  missing  was  the  previous  day’s 
average  substituted.  With  the  data  set  filled  in,  the 
analysis  could  begin.  The  data  files  along  with  the 
BMDP  instruction  files  were  easily  transfered  back  and 
forth  between  the  SSC  and  the  Cyber. 

Discarded  Data 

One  final  point  should  be  made.  As  mentioned 
earlier  in  this  chapter,  only  about  the  last  three  years 
of  the  data  from  the  LANL  tape  were  used.  This  was  due 
to  the  fact  that  the  LANL  tape  contained  the  majority  of 
its  missing  values  (17  out  of  32)  in  the  first  year  of 
the  data  (eg.,  in  1982).  Thus,  to  reduce  the  effect  of 
substitution  of  estimated  missing  values,  this  data  was 
discarded.  In  addition,  since  the  LANL  data  ended  on 
May  6,  1986,  it  was  decided  to  discard  all  data  up  to 
May  8,  1983  which  allowed  the  analysis  of  precisely  the 
last  three  years  of  the  data  (1095  days).  It  was  felt 
that  the  analysis  of  the  most  recent  data  would  allow 
the  development  of  the  most  accurate  and  up  to  date  time 
series  models . 
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In  the  case  of  the  Omni  Tape,  only  the  last  two 


years  were  used  for  the  transfer  function  modeling  (May 


1983  to  April  1985)  because  this  gave  the  maximum 


overlap  of  the  LANL  data.  The  reader  will  recall  from 


earlier  discussions  that  time  synchronized  data  from  two 


series  is  necessary  in  order  to  do  transfer  function 


modeling.  The  overlap  occurred  between  the  dates 


mentioned  above  and  was  more  than  adequate  to  perform 


the  analysis. 


With  regard  to  the  actual  data  which  was  used, 


Appendix  C  is  a  FORTRAN  program  which  was  composed  to 


write  out  all  the  pertinent  data  utilized  in  this  study. 


Appendix  D  is  the  result  of  this  program:  a  complete 


listing  of  the  data. 


Software  and  Computers 


BMDP2T ,  the  "Box-Jenkins  Time  Series  Analysis", 


was  the  program  used  to  perform  all  the  analyses  (Dixon 


and  others,  1985:639-660).  The  initial  step  in  running 


BMDP2T  is  to  compose  a  small  computer  program  containing 


the  necessary  instructions  to  perform  a  basic 


examination  of  the  data  (Dixon  and  others,  1985:640). 


The  stored  data  files  are  accessed  by  either  including 


them  directly  as  part  of  the  instruction  file  (as  called 


for  on  the  Cyber)  or  by  accessing  them  external  to  the 


instruction  file  via  a  data  file  i dent i f i cat i  i 


statement  (as  called  for  on  the  SSC ) .  The  basic 


examination  exhibits  the  data  for  each  time  series  in  a 


graph  of  the  daily  flux  values  over  time  (a  time  plot). 


In  addition,  it  performs  an  initial  calculation  of  the 


autocorrelation  and  partial  autocorrelation  functions  of 


the  data  so  as  to  help  define  an  initial  guess  at  a 


prediction  model.  Such  terms  as  "autocorrelation 


function"  along  with  the  details  of  the  Box  and  Jenkins 


time  series  analysis  methods  will  be  given  later  on  in 


the  Detailed  Methodology  chapter  of  this  paper.  For 


now,  it  is  sufficient  for  the  reader  to  understand  that 


this  type  of  analysis  combines  a  certain  amount  of 


precision  with  an  equal  amount  of  "art"  and  "gut 


feeling".  Said  more  formally,  although  this  method  does 


have  some  "hard  rules",  its  application  to  real  data 


involves  some  interpretation.  Time  series  analysis  is 


an  iterative  and  time  consuming  method.  The  prediction 


models  derived  from  a  time  series  analysis  are  not 


obtained  by  a  definite  set  of  procedures.  Moreover, 


there  may  be  more  than  one  model  which  is  appropriate 


for  the  data  depending  on  the  total  history  or  extent  of 


the  data  which  is  used.  Time  series  models  are  always 


subject  to  updating  when  new  data  becomes  available  to 


add  to  the  history  of  the  series. 
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Two  computer  systems  were  used  in  this  research. 
The  CDC  Cyber  6000  provides  the  fastest  processing 
capability  regardless  of  the  task  (usually  no  more  than 
30  seconds),  but  it  is  a  bit  more  restrictive  than  the 
SSC  Unix  system  which  has  better  manipulative  and  naming 
capabilities  for  files.  The  BMDP  Statistical  Software 
package  is  readily  available  on  both  systems.  The 
current  BMDP  manual  is  the  1985  reprinting  (Dixon  and 
others,  1985).  Documentation  for  the  BMDP  programs  as 
implemented  on  either  system  is  available.  For  the  SSC, 
a  file  named  bmdp . doc  (documentation  file)  is  made 
available  to  the  user  when  the  entire  set  of  BMDP 
programs  is  initially  accessed.  This  file  contains 
helpful  information  concerning  how  to  actually  run  a 
given  BMDP  program.  Included  are  tips  to  increase 
memory  size  (RAM)  slotted  for  each  program  so  that  all 
appropriate  processing  may  be  performed  as  necessary. 
This  one  particular  feature  was  utilized  often  due  to 
the  size  of  the  data  sets  involved  along  with  the  amount 
of  processing.  A  file  similar  to  the  bmdp.doc  file  is 
also  available  for  the  Cyber  implementation  of  BMDP. 

The  on-duty  computer  consultant  can  obtain  a  hard  copy 
of  this  file  for  anyone  desiring  it  as  well  as  show  the 
user  how  to  access  the  BMDP  software  on  either  the  Cyber 
or  the  SSC. 


I V .  Detail ed  Methodol ogy 

The  Box  and  Jenkins  Method 

The  book  entitled  Time  Series  Analysis : 
Forecasting  and  Control  by  G.E.P.  Box  and  G.M.  Jenkins 
(1976)  is  the  backbone  for  the  analysis  in  this  report. 
The  basic  ideas  in  a  time  series  analysis  are  presented 
in  the  first  chapter  of  this  book  (Box  and  Jenkins, 
1976:1-19).  As  compared  to  its  more  classical 
statistical  counterparts,  this  method  is  relatively 
unknown.  Therefore,  a  few  of  the  fundamental  ideas  in 
time  series  analysis  will  be  discussed  along  with  some 
of  the  models  which  form  a  framework  for  a  time  series 
forecast  function.  In  addition,  the  reader  will  be  made 
aware  of  how  a  time  series  transfer  function  is  derived 
Those  with  a  deeper  interest  in  the  workings  of  the  Box 
and  Jenkins  method  are  strongly  urged  to  obtain  a  copy 
of  the  text  for  their  own  perusal. 

Time  Series 

A  series  of  recorded  values  of  a  given  random 
variable  taken  at  equally  spaced  intervals  of  time  is  a 
discrete  time  series .  Since  many  entities  in  the  world 
can  be  termed  random  variables,  analysts  thus  have  the 


opportunity  to  record  periodic  values  of  these  entities 


and  study  them  as  they  behave  over  time  (ie.,  in  a  time- 
series)  . 


A  time  series  analysis  has  three  important 
"areas  of  application"  as  mentioned  in  the  introductory 
chapter.  Two  of  these  have  direct  bearing  upon  this 
study  and  are  noted  as 

(1)  The  forecasting  of  future  values  of  time 
series  from  current  and  past  values. 

(2)  The  determination  of  the  transfer  function 
of  a  system--the  determination  of  a  dynamic 
input-output  model  that  can  show  the  effect  on 
the  output  of  a  system  subject  to  inertia,  of 
any  given  series  of  inputs  (Box  and  Jenkins, 

1976:1) . 

Recalling  the  introduction,  the  reader  will  note  that 
the  objectives  of  this  study  were  to  derive  prediction 
(forecasting)  models  for  energetic  electron  flux  and 
also  to  derive  relationships  (transfer  functions) 
between  the  IMF  Bz»  the  solar  wind  speed,  and  the 
energetic  electron  flux.  This  second  objective  is 
undertaken  with  the  IMF  or  the  solar  wind  speed  viewed 
as  the  input  series  and  the  energetic  electron  flux  in 
any  given  channel  viewed  as  the  output.  As  noted  by  Bo 
and  Jenkins,  a  good  forecast  function  can  provide  the 
basis  for  correct  planning  and  control  (Box  and  Jenkins 
1976:1).  In  this  particular  case,  a  good  forecast 
function  would  allow-  optimal  control  of  spacecraft  and 
sensors  which  may  be  exposed  to  predicted  dangerous 
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levels  of  energetic  electron  flux. 


Also ,  a  valid 
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transfer  function  would  help  to  establish  more  firmly 


the  theoretical  relationship  between  the  IMF,  the  solai 


wind,  and  the  energetic  electron  flux. 


Fundamental s 


A  random  variable  (rv)  of  interest  whose  values 


have  been  recorded  over  equispaced  intervals  of  time 


might  be  designated  as  z .  If  we  denote  the  present  time 


as  t,  then  zt  represents  the  currently  recorded  value 


for  the  rv  of  interest.  Similarly,  one  might  denote  the 


last  recorded  value  of  z  prior  to  the  current  value  as 


zt-i.  The  value  recorded  two  time  periods  prior  to  the 


current  would  be  zt-2  and  so  on.  The  elapsed  time 


period  between  which  the  readings  of  the  variable  are 


made  can  take  on  any  value.  Often  times,  it  will  be 


dictated  by  the  apparatus  used  to  measure  the  variable 


(in  this  case,  the  time  between  satellite  sensor 


measurements  of  the  energetic  electron  flux,  the  IMF 


component  values,  and  the  solar  wind  speed).  The 


important  point  is  that  the  individual  discrete  readings 


of  the  variable  must  be  equispaced  over  time. 


What  is  sought,  then,  is  a  forecast  of  the  value 


of  the  variable  at  some  future  time  t+1.  One  might 


denote  this  forecast  as  £t(l)  where  §  (pronounced 
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zee-hat)  is  the  estimate  of  the  series’  value  at  time 


t+1.  According  to  Box  and  Jenkins,  the  objective  is  to 


derive  a  forecast  function  such  that  "the  mean  square  of 


the  deviations  zt< i  -  5t(l)  between  the  actual  and 


forecasted  values  is  as  small  as  possible  for  each  lead 


time  1"  (Box  and  Jenkins,  1976:2). 


In  addition,  it  is  necessary  to  specify  the 


accuracy  of  each  forecast  function  so  that  "the  risks 


associated  with  decisions  based  upon  forecasts  may  be 


calculated"  (Box  and  Jenkins,  1976:2).  These  accuracies 


are  typically  specified  by  calculating  probability 


limits  (confidence  intervals)  on  either  side  of  each 


forecast  value  (Box  and  Jenkins,  1976:2).  Thus,  if  one 


desired  an  estimate  for  the  energetic  electron  flux  two 


days  in  the  future  based  upon  the  present  history  of  the 


data,  the  forecast  model  employed  might  yield  a  value 


(with  arbitrary  units)  of,  say,  3.6  +  /-  1.2  where  the 


1.2  on  either  side  of  the  estimated  3.6  might  represent 


the  standard  error  for  the  estimate.  This  being  the 


case,  the  given  estimate  along  with  the  stated  +/-  1.2 


standard  error  interval  would  represent  a  confidence 


interval  (or  probability  limits)  for  the  estimate  of 


approximately  68%.  A  95%  confidence  interval  estimate 


would  then  be  represented  by  (approximately)  3.6  +/- 


2.4.  Obviously,  if  one  desires  a  greater  accuracy ,  then 


the  confidence  interval  for  the  estimate  has  to  become 
broader.  In  addition,  the  farther  into  the  future  that 
a  prediction  is  desired,  the  broader  the  confidence 
interval  for  such  an  estimate  will  be. 

To  draw  a  parallel  between  what  has  been 
discussed  thus  far  and  the  data  presented  here  for 
analysis,  Figure  1  is  shown  below.  In  it  is  a  plot 
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Fig  1.  Hypothetical  Time  Series  Plot 
(Box  and  Jenkins,  1976:2) 


of  a  hypothetical  time  senes  for  energetic  electron 
flux  readings  recorded  over  a  period  of  30  days.  If  the 
present  time  is  t,  then  what  is  desired  is  an  estimate 
of  the  series'  value  at  time  t+1  where  1  is  perhaps  two 
days  hence.  Note  that  the  confidence  intervals  shown 
illustrate  that  the  farther  into  the  future  one 
estimates,  the  broader  the  confidence  interval  becomes. 
Likewise,  the  greater  the  specified  accuracy  of  the 
estimate  at  any  desired  time  in  the  future,  the  broader 
the  interval . 

Notice  that  the  word  "probability"  is  used  in 
the  above  discussion.  This  is  because  a  time  series 
model  is  a  probabilistic  or  stochastic  model  as  opposed 
to  a  deterministic  model.  In  a  stochastic  model,  one 
calculates  "the  probabi I i ty  of  a  future  value  lying 
between  two  specified  limits"  (Box  and  Jenkins,  1976:7). 
It  is  also  important  to  note  that  the  actual  time  series 
recorded  is  just  one  of  an  infinite  number  of  possible 
series  which  could  have  been  generated  by  the  given 
stochastic  model  for  the  series  (Box  and  Jenkins, 

1976:7) . 


Categories  of  Processes 


Generally,  there  are  two  different  classes  of 
time  series  models.  Stationary  models  derive  their  name 


from  the  fact  that  the  series  they  represent  tend  to  be 


in  equilibrium  about  a  constant,  never-changing,  mean 
level  (Box  and  Jenkins,  1  976:7).  A  nonstat i oriary  model, 
however,  is  used  to  represent  a  series  (also  called  a 
process)  where  there  is  no  natural  mean  or  equilibrium 
point  (Box  and  Jenkins,  1976:7).  As  will  be  explained 
later,  a  given  time  series  usually  reveals  whether  or 
not  it  is  stationary  when  its  autocorrelation  function 
is  plotted.  If  the  autocorrelations  decrease  quickly  as 
the  lag,  k,  increases,  the  series  is  considered 
stationary.  In  all  the  cases  studied  here,  the  series 
were  stationary.  A  particularly  wide  class  of  processes 
which  "provides  a  range  of  models,  stationary  and 
non-stationary ,  that  adequately  represent  many  of  the 
time  series  met  in  practice"  is  known  as  the 
autoregressive-integrated  moving  average  processes  or 
ARIMA  processes  (Box  and  Jenkins,  1976:8). 


Operators  and  Assumptions 

In  describing  the  primary  time  series  models,  it 
is  necessary  to  discuss  a  couple  of  the  basic  operators 
used.  Perhaps  the  two  most  important  are  the  backward 
shift  operator  denoted  by  B  and  the  backward  difference 
operator  denoted  by  V  (pronounced  "del").  These  two 
operators  are  defined  as  follows: 
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Bz  »  =  z  t -  1 


Bmz  t  =  z  t - b 


Vz %  -  z  t  -  z  t  -  1 


=  (1  -  B ) z 


VmZ  X  =  zt  -  Zt-1 


Note  that  the  backward  difference  operator  can  be 


represented  in  terms  of  the  backward  shift  operator. 


Also,  the  parameter  m  may  be  thought  of  as  the  "power" 


of  the  operator. 


In  addition  to  these  basic  operators,  the  time 


series  models  used  by  Box  and  Jenkins  view  the  output  of 


any  given  time  series  as  being  generated  from  an  input 


series  of  independent  shocks,  at,  with  the  shocks  being 


random  drawings  from  a  fixed  normal  distribution  of  mean 


zero  and  variance  <7a  2  .  "Such  a  sequence  of  random 


variables  at,  at-i,  at-j, . . .  is  called  a  white  noise 


process  by  engineers"  (Box  and  Jenkins,  1976:8).  When 


the  model  for  the  series  under  study  is  finally  fitted, 


the  white  noise  represents  the  res  idual s  or  the 


differences  between  the  actual  series  values  and  those 


which  the  model  would  predict  at  the  same  points  in  time 


at  which  the  actual  series  value's  were  recorded.  This 
sequence  of  random  white  noise  numbers  is  transformed 
into  the  observed  output  time  series,  zt,  by  another 
operator  known  as  a  linear  filter,  0(B)  where  B  is  again 
the  backshift  operator  and  the  notation  0(B)  indicates 
that  varying  powers  of  backshifting  are  used  throughout 
the  operator.  For  example,  one  might  represent  a  time 
series  of  interest  by  the  model 

zt  =  M  +  0(B)at  ( 1 ) 


where  the  filter,  0 ( B ) ,  is 

0(B)  =  1  +  0,B  +  02B*  + 


A  finite  or  infinite/convergent  set  of  parameters  0i  , 
02,  ...  is  said  to  classify  the  process  as  stationary. 

In  this  case,  the  value  [1  is  thus  the  natural  mean  or 
equilibrium  point  for  the  process  under  study.  If  the 
set  of  parameters  is  infinite  and  not  convergent 
however,  then  the  process  is  said  to  be  nonstationary, 
and  in  this  case,  jl  is  only  a  "reference  point  for  the 
level  of  the  process"  (Box  and  Jenkins,  1976:8-9). 


Spec  i fie  Models 


An  autoregress ive  mode  1  for  a  time  series  is 
expressed  "as  a  finite,  linear  aggregate  of  previous 
values  of  the  process  and  a  shock  at"  (Box  and  Jenkins, 
1976:9).  If  the  values  of  a  process  are  zt,  Zt-i,  zt-2, 

. . . .  as  before  and  the  deviations  of  each  of  these 
values  from  the  mean,  /J,  are  represented  by  5ft,  zt-i, 
zt-2,  ....  then  an  autoregressive  (AR)  model  of  order  p 

is  represented  by 

2ft  =  0i  5f  t  -  i  +  02  z  t  -  2  +  ...  +  0p  zt-p  +  at  (2) 

which  as  the  reader  will  note,  is  simply  a  regression 
equation  where  the  value  5f  t  is  regressed  on  previous 
values  of  itself.  The  term  at  is  analogous  to  the  error 
term  in  a  regression  and  in  this  context  represents  the 
random  shock  input  in  the  current  time  period.  Thus, 
the  reason  for  the  name  autoregressive  becomes  apparent 
(Box  and  Jenkins,  1976:9).  The  autoregressive  operator 
may  then  be  defined  as 

0(B)  =  1  -  0i  B  -  02B*  -  ...  -  0PB p 

and  in  doing  so,  the  autoregressive  model  above  may  be 
rearranged  and  "written  economically"  so  that 

0(B)zt  =  at  (Box  and  Jenkins,  1976:9) 


42 


As  Pox  and  Jenkins  point  out.,  it  is  not  difficult  to  see 
that  the  autoregressive  model,  (2),  is  a  special  ease  of 
the  linear  filter  model,  (1),  discussed  previously  (Box 
and  Jenkins,  1976:10).  To  show  this,  the  z  terms  on  the 
right,  hand  side  of  (2)  may  be  progressively  replaced 
with  their  expanded  values  according  to  the  model.  For 
instance,  the  term  ?t-i  on  the  right  side  of  (2)  may  be 
replaced  and  thereby  eliminated  by  substituting 


Zt-l  r  0lZt-2  +  <t>2Z  t  -  3  +  ...  +  0p? t - p - i  + 


Similar  substitutions  for  Zt-2,  and  so  on  will 
eventually  yield  an  infinite  series  in  the  a’s  (Box  and 
Jenkins,  1976:10).  Having  made  all  the  substitutions, 
one  can  thus  see  that  not  only  is 


0(B)  z  t  =  at 


as  already  pointed  out,  but  also 


z  t  =  0  (  B  )  a  t 


Therefore,  it  is  obvious  that 


0(B)  =  0_1(B)  (Box  and  Jenkins,  1976:10) 


The  reader  will  note  that  the  utility  of  using  0(B)  over 
0(B)  in  model  development  stems  from  the  fact  that  while 


rAVwv, 


0(B)a 


For  discrete  data  such  as  this  study  will 
analyze,  the  trans f er  f unction  between  an  output  time 
series,  Yt,  and  an  input  time  series,  Xt|  may  be 
effectively  represented  by  a  difference  equation 
involving  backward  shift  operators  which  operate  on  both 
the  input  and  the  output  series  (Box  and  Jenkins, 
1976:13).  This  type  of  representation  is  given  as 


(  1  -  6iB  <5  r  B  r  )  Y  t 

=  (  Do  -  D  i  B  tdsB  s  )  X  t  -  b  +  a 


=  (  (J0Bb  -  6JiBb  ♦  1  -  .  .  . 


-  lUsBb  ♦  s)Xt  +  at  (  8  ) 


<5  ( B  )  Y  t  =  0>(B}B»Xt 

=  0(B)X» 


where  the  parameter  b  is  the  number  of  time  periods 
(dead-time  or  pure-delay)  between  the  input  value  and 
its  observed  effective  output.  The  parameters  r  and  s 
express  the  order  of  the  left  side  and  right  side 
operators  in  the  above  difference  equation  and  at  is  a 
random  noise  component.  Another  way  of  stating  the 


relationship  is  that  Y*  and  Xt  are  linked  by  a  linear 
filter  such  that 

Y  t  =  \'oX  t  +  v  i  X  t  -  1  +  v  jX  i  -  2  +  ...  +  at  (  9  ) 


or 


Yt  =  v  ( B  )  X  »  +  at 

where  v(B)  is  the  transfer  function  and  can  be  stated 
more  explicitly  as 

v ( B )  =  vo  +  v  i  B  +  v2B2  +  ...  (10) 

Thus  v ( B )  is  simply  a  ratio  of  the  right  side  operator 

of  order  s  and  the  left  side  operator  of  order  r  in  the 
difference  equation  (8).  The  reader  should  take  note  of 
the  fact  that  for  an  input-output  relationship  in  which 
the  dead-time  between  the  input  and  its  observed  output 
is  equal  to  b  time  periods,  the  first  b  weights  in 
equation  (10)  (eg.,  Vo,  vj  ,  ...  ,  vt>-i)  are  zero  (Box 

and  Jenkins,  1976:14). 

With  regard  to  the  data  available  for  this 
study,  it  is  hoped  that  a  model  of  the  types  discussed 

here  will  fit.  The  tasks  thus  implied  with  respect  to 

th€'  objectives  listed  in  the  introduction  to  this  report 


are  as  foil ows : 


(1)  Develop,  if  possible,  a  time  series  model  of  the 
type  AR ,  MA ,  or  ARMA  for  each  of  the  energetic  electron 
channels  discussed. 

(2)  Use  the  developed  models  to  make  forecast. s  of  the 
flu\  values  and  compare  these  forecast  values  to  the 
actual  values  to  assess  the  model’s  utility  in 

predi ction. 

(3)  Develop  models  of  the  same  type  for  solar  wind  and 
IMF  Bz  in  an  attempt  to  show  that  the  solar  wind  speed 
or  the  IMF  Bz  (considered  as  inputs)  are  related  to  the 
energetic  electron  flux  levels  (considered  as  outputs) 
via  a  transfer  function. 


Deduc i ng  the  Mode  1 

The  "precise  steps"  to  be  followed  in  deriving 
these  models  have  so  far  not  been  mentioned.  As  stated 
in  the  Chapter  III,  the  Box  and  Jenkins  method  combines 
a  certain  amount  of  precision  with  an  equal  amount  of 
art.  The  "precision"  in  deriving  the  appropriate 
model (s)  comes  with  the  determination  of  the 
autocorrelation  and  partial  autocorrelation  functions . 
The  "art"  comes  in  the  form  of  how  well  one  can  deduce 
the  type  of  model  from  these  special  functions. 

Deducing  the  model  is  an  iterative  and  time  consuming 
process  which  involves  close  scrutiny  of  these  functions 
as  many  models  are  fitted  to  the  data  until  eventually 
(and  hopefully)  one  of  the  fitted  models  leads  to  a 


reasonably  small  set  of  residuals. 

Thus,  a  necessary  first  step  is  the 
determination  of  the  autocorrelation  functions  for  the 


r^n' 


> 


t 


time  series  of  interest.  For  a  series  with  N  discrete 
values,  the  kth  lag  autocorrelation  is  defined  as 

r  k  =  Ck/Co  (11) 


where 

1  N  -  k 

Ck  =  -  L(zx  -  z)(Zfk  -  z),  k  =  0,1,2,..  .,K  (12) 

NT  t  =  i 

is  an  estimate  of  the  autocovariance  at  lag  k,  and  zT  is 
the  mean  of  the  time  series  (Box  and  Jenkins,  1976:32). 
The  lag  (k)  has  to  do  with  the  number  of  time  periods 
between  recorded  values  in  the  process.  For  instance, 
if  the  series  of  interest  consists  of  1000  daily 
readings  (eg.,  readings  taken  precisely  24  hours  apart), 
then  the  autocorrelation  at  lag  k  =  2  days  would  consist 
of  the  sum  of  all  autocovariances  of  readings  which  lag 
each  other  by  exactly  two  days  (C2)  divided  by  the  sum 
of  all  autocovariances  which  lag  each  other  by  exactly 


I 


zero  days  (Co).  The  reader  may  also  note  that  Co  is 
nothing  more  than  the  simple  variance  of  all  the 
observations.  Thus,  the  autocorrelations  for  any  lag 
from  zero  on  up  to  N  -  1  may  be  calculated  by  the 
formulas  (11)  and  (12).  The  display  of  these 
autocorrelations  in  a  bar  chart  showing  the  computed 
values  plotted  against  progressively  higher  lags  forms 
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what  is  commonly  refered  to  as  the  autocorrelation 
f  unct i on . 

Figure  2  is  an  example  of  such  a  chart.  Tin- 
autocorrelation  function  along  with  a  nearly  identical 
illustration  of  the  partial  autocorrelation  function 
(Figure  3),  form  the  fundamental  tools  in  the  iterative 
procedure  of  identifying  the  model  which  fits  a  given 
time  series.  Whereas  the  shape  of  the  autocorrelation 
function  allows  a  guess  as  to  the  type  of  model  which 
should  be  entertained  (eg.,  AR,  MA,  ARMA,  etc.)  the 
parti a 1  autocorrelation  f unc t i on  is  used  to  help 
determine  the  order  of  the  model.  Box  and  Jenkins  liken 
the  use  of  the  partial  autocorrelation  function  to  that 
of  "deciding  on  the  number  of  independent  variables  to 
be  included  in  a  multiple  regression"  (Box  and  Jenkins, 
1976:64).  They  also  include  a  detailed  discussion  of 
how  to  obtain  the  partial  autocorrelation  function  (Box 
and  Jenkins,  1976:64-65).  Basically,  the  partial 
autocorrelation  function  shows  the  autocorrelation  at 
higher  lags  after  the  effects  of  variables  up  to  a 
certain  lag  have  been  regressed  out.  Therefore,  the 
partial  at  lag  4  is  the  autocorrelation  at  lag  4  having 


accounted  for  the  effects  of  lag  1,  2,  and  3. 


Fig  2.  Example  Autocorrelation  Function 


An  example  to  clarify  the  above  discussion  is 
justified.  Referring  to  Figure  2,  let  the  assumption  be 
that  the  autocorrelation  function  shown  represents  a 
time  series  which  is  under  study.  One  can  see  that  the 
locus  of  points  formed  by  an  imaginary  curve  connecting 
the  tops  of  the  charted  values  has  a  shape  which  is 
characteristic  of  some  sort  of  exponential  function. 
According  to  Box  and  Jenkins,  this  is  a  classic 


One  point  should  be  made  before  determining  the 
order.  Box  and  Jenkins  give  formulas  for  calculating 
the  standard  error  of  the  autocorrelation  or  partial 
autocorrelation  at  any  chosen  lag  (Box  and  Jenkins, 
1976:34,65).  Twice  the  standard  error  is  chosen  as  the 
level  of  significance  for  any  calculated  autocorrelation 
or  partial  autocorrelation  value.  In  other  words,  if 
the  computed  correlation  value  for  any  given  lag  k  does 
not  exceed  twice  the  standard  error,  then  that 
correlation  is  effectively  regarded  as  having  a  value  of 
zero.  Again  referring  to  Figure  2  and  3,  the  dotted 
lines  on  each  bar  chart  reflect  twice  the  standard  error 
at  each  lag.  Observe  that  even  though  noticeable  values 
have  been  computed  past  the  second  lag  on  the  partial 
autocorrelation  plot,  in  terms  of  the  indicated  standard 
error  they  are  insignificant  and  thus  are  regarded  as 
zero.  In  viewing  the  partial  autocorrelation  function 
(Figure  3),  one  can  see  that  the  function  shows  a 
definite  cutoff  after  the  lag  k  =  2  indicating  that  the 
process  being  analyzed  is  very  likely  a  second  order 
series.  This  being  the  case,  the  model  would  be 
initially  estimated  as  AR ( 2 )  and  would  thus  have  the 


Similar  procedures  are  followed  for  all  lime  series. 


The  shape  and  character  of  the  autocorrelation  and 
partial  autocorrelation  functions  are  thus  the 
determining  factors  in  the  initial  guess  as  to  the  model 
type.  Box  and  Jenkins  devote  chapter  6  to  the 
discussion  of  how  to  identify  the  type  of  model  along 
with  estimating  its  parameters  (0i  and  02  in  this  case) 


(Box  and  Jenkins,  1976:171-207). 

The  reader  may  now  be  wondering  how  the 
parameters  of  the  entertained  model  are  calculated.  As 
already  stated,  a  significant  amount  of  detail  is 
contained  in  the  text  regarding  model  parameter 
estimation  (Box  and  Jenkins,  1976:187-193,208-284).  Let 
it  suffice  to  say  that  a  modern  statistical  software 
package  such  as  BMDP  allows  the  estimation  of  these 
parameters  with  a  minimum  of  effort  (Dixon  and  others, 


1985:644-645).  The  main  task  of  the  analyst  is  to 
merely  decide  on  the  type  and  order  of  model  he/she 


would  like  to  fit  to  the  data. 


Model  Checkini 


After  the  model  is  identified  and  its  parameters 


estimated,  the  model  must  then  be  tested  to  see  whether 
a  good  fit  has  indeed  occurred.  If  the  fit  is  good, 


then  the  model  may  be  retained  for  forecasting.  If  not, 
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then  a  new  model  must  be  entertained  and  likewise 
tested.  This  process  continues  until  a  good  fit  is 
obtained.  Box  and  Jenkins  refer  to  the  testing  of  the 
model  as  "model  diagnostic  checking"  (Box  and  Jenkins, 
1976:285-299).  One  of  two  methods  they  describe  to  do 
diagnostics  is  that  of  displaying  "the  autocorrelation 
function  of  the  residuals"  (Box  and  Jenkins,  1976:285). 

Earlier  in  this  chapter,  it  was  stated  that  the 
time  series  models  of  Box  and  Jenkins  employ  a  white 
noise  process  (a  series  of  normal  shocks  with  mean  zero 
and  variance  <Ja2).  In  so  doing,  this  white  noise  is 
transformed  into  a  model  which  adequately  represents  the 
series  under  study  by  linear  filtering.  The  model  thus 
derived  can  be  rearranged  so  that  estimated  values  of 
the  originally  input  white  noise  process  maj'  then  be 
computed.  When  this  is  done,  the  estimated  white  noise 
(or  shocks)  are  designated  as  at  and  are  also  called  the 
residuals  (Box  and  Jenkins,  1976:289).  The  plot  of  the 
autocorrelation  function  of  these  residuals  is  one  of 
the  primary  methods  used  to  diagnose  the  adequacy  of  the 
model.  If  the  autocorrelations  of  the  residuals  are 
significant  at  any  lag,  then  a  better  model  can  most 
likely  be  found.  The  text  suggests  ways  to  use  the 
residual  autocorrelations  to  implement  a  better  model 


(Box  and  Jenkins,  1976:298-299). 


If  they  are 


insignificant,  then  the  model  is  adequate  and  may  be 
used  for  forecasting.  Again,  twice  the  standard  error 
determines  the  significance  of  the  computed 
autocorrelation  at  each  lag.  This  method  is  easily 
accomplished  with  the  BMDP  package  and  used  exclusively 
in  this  study. 

To  recap,  the  building  of  a  time  series  model  is 
an  iterative  process  (Box  and  Jenkins,  1976:19).  In 
most  cases,  a  five  step  process  takes  place: 

(1)  The  general  class  of  models  is  postulated.  A  plot 
of  the  actual  time  series  may  reveal  the  type  of  model 
to  the  "practiced  eye". 

(2)  The  initial  model  to  be  entertained  is  then 
identified  from  the  autocorrelation  and  partial 
autocorrelation  functions. 

(3)  The  parameters  of  the  entertained  model  are 
estimated . 

(4)  The  model  is  diagnostically  checked  by  examining 
the  autocorrelation  of  the  residuals  computed  from  the 
newly  fitted  model  with  its  estimated  parameters. 

(5)  If  the  model’s  residuals  are  insignificant,  the 
model  may  be  used  for  forecasting.  If  not,  a  new  model 
must  be  suggested  based  on  the  residual  autocorrelation 
plot  examination  and  steps  2  through  5  reaccomplished 
until  a  suitable  model  is  found. 


Transfer  Funct i on  Mode 1 s 

As  discussed  earlier,  a  transfer  function  links 
the  output  time  series  of  a  given  dynamic  system  with 
the  input  time  series  to  that  system.  From  this 
perspective,  one  may  be  better  able  to  understand  the 
actual  dynamics  of  the  system  in  a  purely  mathematical 
sense  if  not  in  the  physical  sense.  Such  mathematical 


k.  .  , 
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insights  are  almost  certain  to  eventually  lead  to  bettei 
physical  understanding  as  well. 

Box  and  Jenkins  go  into  considerable  detail  on 
the  subject  of  transfer  function  modeling  in  chapters  10 
and  11  of  their  text  (Box  and  Jenkins,  1976:335-420). 

In  pei haps  their  best  illustration  of  how  a  transfer 
function  model  is  formed,  they  exhibit  the  necessary 
steps  in  an  example  problem  concerning  the  output  of 
carbon  dioxide  (CO2)  from  a  gas  furnace  which  has 
influential  inputs  of  air  and  methane  ( CH 4  )  (Box  and 
Jenkins,  1976:381-386).  In  a  very  fortunate 
coincidence,  the  BMDP  manual,  in  discussing  their  own 
example  problem  to  illustrate  transfer  function 
miodeling,  uses  the  same  gas  furnace  problem  of  Box  and 
Jenkins  (Dixon  and  others,  1985:650-655).  Thus  the 
reader  is  able  to  match  the  logic  used  in  forming  the 
computer  solution  to  the  problem  with  the  theoretical 
logic  contained  in  Box  and  Jenkins.  For  this  reason, 
further  discussion  of  transfer  function  modeling  will  be 
limited  to  the  main  tool  of  analysis  in  this  method 
which  is  the  cross  correlation  f uncti on .  The  reader 
desiring  more  detail  on  transfer  function  modeling  is 
therefore  urged  to  obtain  a  copy  of  Box  and  Jenkins  and 


the  BMDP  manual. 


V4** 


Just  as  the  autocorrelation  and  partial 


autocorrelation  functions  yield  much  information  about  a 


given  single  time  series,  the  cross  correlat  i  cj  r>  f  uriet  i  on 


reveals  much  about  the  relation  between  anv  two  time 


series.  A  cross  correlation  between  two  time 


series  may  be  appropriate  at  any  time  one  believes  that 


the  two  interact  in  some  sort  of  a  dynamic  system.  As 


brought  out  in  the  literature  review,  this  is  the  belief 


concerning  the  energetic  electron  flux  levels  and  the 


solar  wind  speed/IMF  B2.  Formulae  for  determining  the 


cross  correlation  at  any  desired  lag  between  any  two 


dynamically  related  time  series  are  contained  in  the 


text,  and  the  reader  will  notice  that  they  are  very 


similar  to  those  given  for  determining  the 


autocorrelations  (Box  and  Jenkins,  1976:374).  These 


formulae  are  as  follows: 


r  x  y ( k  )  -  Cxy(k)/SxSy 


where 


c  x  >  (  k  )  = 


1  n  -  k 

-  £(x t  -  X) (y t ♦ k  -  y )  k  =  0,1,2,...  (14) 


n  »  i 


c  X  >■  (  k  )  = 


1  n  ♦  k 

X  ( y  t 


y)  (x,  -k  -  X)  k  =  0,-1 , -2, 


In  equaLion  (13),  i  ,,(k)  is  the  cross  roi-rrlalion  al  lag 
k  between  an  input  series  x  and  an  output  series  y.  The 
term  c i >  is  the  cross  covariance  between  the  two  scries, 
and  equations  for  obtaining  it  are  shown  in  (14)  and 
(15)  where  n  stands  for  the  number  of  pairs  of  series 
values  ( x  t | y  t )  and  x  and  “  are  the  respective  series 
means . 

Once  the  cross  correlations  at  each  lag  are 
determined,  they  are  plotted  in  a  bar  chart  which  is 
very  similar  to  Figures  2  and  3  shown  earlier  in  this 
chapter.  This  plot  is  referred  to  as  the  cross 
correlation  function .  From  this  chart  the  model  is 
tentatively  identified.  An  estimation  of  the  dead-time 
parameter,  b,  may  be  made  directly  from  the  crosa 
correlation  function  along  with  estimates  of  the 
parameters  r  and  s  from  the  difference  equation  (8). 

The  parameter  b  is  usually  set  equal  to  the  value  for 
the  first  positive  lag  at  which  a  significant  cross 
correlation  occurs.  Thus,  if  the  first  significant 
cross  correlation  for  two  dynamically  connected  series 
occurs  at  lag  k=2,  then  b  is  set  equal  to  2.  The  model 
estimation  and  final  identification  process  then 
proceeds  iteratively  as  it  did  for  the  single  series 


mode  1 s . 


One  final  point  which  should  be  made  is  that  the¬ 


me  t  hod  of  transfer  function  modeling  used  by  BMDP2T  is 
that  of  "prewhitening  the  input"  (Box  and  Jenkins, 
1976:379-380).  This  method  essentially  recognizes  that 
two  series  which  are  dynamically  related  also  have  an 
element  of  noise  involved.  Thus,  any  output  series 
modeled  from  a  related  input  series  is  influenced  not 
only  by  the  input  but  by  noise  as  well.  The  amount  of 
noise  present  in  the  model  can  thus  be  an  indicator  of 
how  deeply  related  any  two  series  are.  Said  another 
way,  the  greater  the  noise  component,  the  less  influence 
the  input  has  on  the  output  series.  This  noise 
component  is  modeled  in  the  transfer  function  models 
determined  in  this  study.  Also,  and  perhaps  more 
importantly,  the  prewhitening  method  is  used  since  it 
can  provide  initial  rough  estimates  of  the  vi,  parameters 
of  the  transfer  function  equation  (10)  (Box  and  Jenkins, 
1976:380).  With  these  initial  rough  estimates,  the 
model  fitting  procedure  can  usually  be  implemented  in  a 
more  efficient  manner. 


BMDF  Examples 

Appendix  E  of  this  report  contains  three  example 
time  series  analysis  problems  as  performed  using  program 


✓ 


BMDP2T.  Example-  1  was  done'  on  the-  SSC,  while?  Example  2 
and  Example  3  were  done  on  the  CYBER. 

Example  1  is  typical  of  the  early  stages  of 
identification  of  a  model  for  any  given  time  series.  In 
this  particular  run,  energetic  electron  channel  SEEI  is 
analyzed,  and  the  plot  of  the  autocorrelation  and 
partial  autocorrelation  functions  indicate  a  stationary 
series  with  autoregressive  tendencies.  In  explanation, 
the  autocorrelation  function  decreases  rather  quickly  as 
the  lag,  k,  increases,  while  the  partial  autocorrelation 
function  has  a  definite  cutoff  after  lag  k  =  2  (page 
138).  This  would  indicate  a  probable  best  fit  of  an  AR 
model  with  first  and  second  order  terms.  However,  the 
analyst,  not  being  seasoned  in  the  "art"  of  time  series 
analysis,  decided  to  attempt  to  fit  an  AR  mode)  with 
first  and  27th  order  terms  as  shown  in  the  ARIMA 
paragraph  on  page  139  of  Appendix  E.  On  the  page  140, 
the  results  of  the  attempted  fit  yield  a  model  with 
estimated  parameters  of  0.7856  and  0.1039  (by  the 
backcasting  method).  The  reader  will  notice  that  the 
T-RATIO  of  the  first  order  is  very  high  (23.16) 
indicating  that  the  model  should  most  definitely  contain 
a  first  order  AR  term.  However,  the  T-RATIO  of  the  27th 
order  is  only  1.91  indicating  that  perhaps  this 
parameter  is  inappropriate.  Nevertheless,  if  one 


believes  in  the  validity  of  the  parameters  as  estimated, 
the  model  would  be 


(1  -  0.7856B  -  0 . 1  039B  2  7  )  z?  t  =  a, 

where  zt  represents  the  current  estimated  deviation  of 
SEEI  from  the  series  mean.  Using  the  plot  of  the 
autocorrelations  of  the  residuals,  at,  resulting  from 
the  attempted  fit  of  the  model  (page  141),  the  reader 
can  see  that  the  fit  is  inappropriate  since  the 
residuals  show  significant  (ie.,  greater  than  two 
standard  errors)  autocorrelations  at  lags  k  =  1  ,  2 ,  and 
4.  Thus,  the  analyst  should  try  to  fit  a  model  which 
better  explains  the  data.  The  best  guess  of  the  author 
would  be,  as  stated  earlier,  an  AR  model  with  first  and 
second  order  terms.  This  iterative  identification 
process  is  an  illustration  of  the  "art"  of  time  series 
analysis.  The  person  with  a  better  of  grasp  of  this 
kind  of  analysis  will  be  able  to  converge  to  a  more 
appropriate  model  in  a  fewer  number  of  tries. 

Example  2  (page  144)  is  much  the  same  as  Example 
1.  Here  the  series  under  study  is  energetic  electron 
channel  SEETII.  The  analyst  has  already  examined  the 
autocorrelation  and  partial  autocorrelation  functions  of 
the  series  and  thus  has  decided  to  spare  computer 
processing  time  by  not  including  the  autocorrelation 
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functions  in  the  analysis. 


In  the  ARIMA  paragraph  on 


page  145,  an  attempt  is  made  to  fit  an  ARMA  model  with 
AR  terms  of  first  and  27  th  order  along  with  an  MA  term 
of  first  order.  Using  the  parameters  estimated  by  the 
backcasting  method,  the  model  is 

(1  -  0.9850B  -  0.1701B27)?,  =  (1  +  0.1227B)a, 

The  fit  of  this  model  is  somewhat  better  than  that  of 
Example  1  as  evidenced  by  the  lesser  number  and 
magnitude  of  the  significant  autocorrelations  of  the 
residuals  (page  147).  The  analyst  may  want  to  accept 
the  fitted  model  or  reiterate  and  attempt  to  fit  a 
different  one.  A  forecast  of  future  values  is  also  made 
(page  148)  based  upon  the  model  as  estimated  above.  The 
reader  can  see  that  the  standard  error  of  the  forecast 
is  0.02037  for  the  initial  day.  It  should  also  be  noted 
that  the  standard  error  increases  as  one  attempts  to 
forecast  farther  into  the  future.  The  analyst  may  also 
use  the  forecast  paragraph  to  check  the  accuracy  of  the 
fitted  model  by  requesting  that  the  forecasting  take 
place  over  a  period  for  which  actual  series  values  are 
available.  In  the  example,  days  345  through  365  have 
actual  values  which  may  be  compared  to  the  forecast 
values. 
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Example  3  is  a  transfer  function  model 
visualizing  the  IMF  Bz  component  as  the  input  series  and 
tiie  SEESSD  channel  as  the  resulting  output.  The 
important  point  to  note  here  is  that  the  plot  of  the 
cross  correlations  (page  153)  shows  only  one  significant 
value  at  lag  k  =  1.  Thus,  the  analyst  might  surmise  a 
very  simple  transfer  function  model  with  the  delay 
parameter,  b,  set  equal  to  1  and  parameters  r  and  s  (the 
left  side  and  right  side  orders  of  the  difference 
equation  (8))  both  equal  to  zero.  Hence,  (r,s,b)  would 
be  (0,0,1).  The  analyst  would  then  use  the  procedure 
described  in  Chapter  10  of  Box  and  Jenkins  to  derive 
initial  estimates  of  of  the  operators  6(B)  and  CJ(B) 

(Box  and  Jenkins,  1976:345-351).  With  these  initial 
estimates,  BMDP2T  is  then  run  interactively  in  a  series 
of  trial  and  error  steps  until  a  suitable  model  is 
identified  (Dixon  and  others ,  1985:650-656).  The 
iterative  part  of  developing  a  transfer  function  model 
is  mainly  contained  in  deriving  a  good  model  for  the 
input  series  and  for  the  noise  component.  Also,  as 
discussed  in  Box  and  Jenkins,  there  may  be  more  than  one 
transfer  function  model.  That  is  to  say,  r,  s,  and  b 
may  be  set  equal  to  other  values  and  other  models  may  be 
fit  with  equally  good  (or  poor)  results  (Box  and 
Jenkins,  1976:387).  The  final  transfer  function  model 


derived  here  and  used  to  obtain  the  forecasts  on  page 
161  is 

Yt  =  -24.0 1BX  t  +  (1/(1  -  0.6125B  -  0 . 1  9 7 3B 2 7 ) ) a t 

The  terms  in  parentheses  represent  the  noise  portion  of 
the  model.  The  first  term  on  the  right  hand  side 
represents  that  portion  of  the  output  series  (where  Yt 
is  the  output  or  SEESSD  in  this  case)  which  may  be 
explained  directly  by  the  input  series  (where  Xt  is  the 
input  series  or  IMF  Bz). 

Again,  the  interested  reader  is  strongly 
advised  to  obtain  a  copy  of  both  the  text  and  the  BMDP 
manual  in  order  to  fully  grasp  the  time  series  analysis 
procedures.  The  results  of  the  Box  and  Jenkins  method 
as  applied  to  the  data  received  will  now  be  presented. 


V .  Resul  t.s 


Ind i v i dua I  Channel  Models 

Tables  I  -  IV  beginning  on  the  next  page 
summarize  the  time  series  analysis  results  for  each  of 
the  energetic  electron  channels.  In  other  words,  the 
series  of  discrete  daily  average  values  from  each 
channel  (SEESSD,  SEEI,  SEEII,  and  SEEIII)  were 
separately  analyzed  in  order  to  try  to  fit  an 
appropriate  model  to  each  individual  channel . 

The  first  thing  the  reader  will  notice  is  that 
no  individual  models  for  SEEIV  are  included.  The  reason 
for  this  is  simply  that  the  data  in  the  SEEIV  channel 
did  not  lend  itself  to  this  type  of  analysis.  SEEIV 
contains  recorded  values  of  energetic  electron  count 
rates  in  the  range  of  9.7  -  16  MeV.  An  initial 
examination  of  the  SEEIV  data  over  the  last  765  days 
(from  April  84  until  May  86)  revealed  that  the 
autocorrelations  were  effectively  zero  at  all  lags  save 
one:  At  lag  k  =  35,  the  computed  autocorrelation  was 

0.199  which  was  significant.  Likewise,  the  computed 
partial  autocorrelations  were  all  effectively  zero 
except  at  lag  k  =  35  where  the  computed  value  was  once 
again  0.199.  The  fact  that  nothing  could  be  said  about 
the  SEEIV  data  due  to  the  lack  of  any  identifiable 


TABLE  I 


NV 


Suggested  Models  for  the  SEESSD  Channel  (1.2  -  1.8  MeV) 


Data  Length 

Suggested  Model 

S  .E  . 

100  days 
(Jan  86-May 

86  ) 

(1  -  .  8250B 
=  (1  + 

-  .  1  7  76B  2  7 
. 4 4 0 IB ) a  t 

)  z  t 

157 .49/ 
174.29 

365  days 
(May  85-May 

86  ) 

(1  -  . 7  4  70B 
=  (1  + 

-  . 2329B  2  7 
.  2 7 2 7B ) a  t 

)  z  t 

227 . 34/ 
241.04 

Approximate  series 

mean  =  224 

.  49 

pattern  in  the  autocorrelations  meant  that  this  data,  by 
itself,  could  not  be  analyzed  via  time  series  analysis. 
However,  transfer  function  modeling  of  this  channel  was 
possible  and  will  be  shown  later.  The  reader  will 
recall  that  the  autocorrel at i on  functions  are  the 
primary  tool  in  identifying  a  series  model.  With  no 
values  other  than  zero,  identification  becomes 
impossible.  Sample  autocorrelation  and  partial 
autocorrelation  functions  for  all  the  energetic  electron 
channels  are  included  in  Appendix  F. 

In  examining  Table  I,  the  reader  will  notice 
that  a  couple  of  different  models  are  identified.  The 
first  model  is  an  AR1,27  MAI  representation  of  the 
SEESSD  channel  over  the  last  100  days  of  the  available 
data  (27  Jan  86-6  May  86).  The  second  model  is  a 
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similar  representation  of  the  last  year  of  the  data  (7 
May  85  -  6  May  86).  These  two  models  represent  the 
"best  fits"  that  the  author  was  able  to  identify  it, 
terms  of  (1)  the  least  residual  mean  square  values  and 
(2)  the  smallest  forecast  standard  errors.  In  both 
cases,  the  plot  of  the  autocorrelation  function  of  the 
residuals  showed  a  good  fit  (no  significant  correlations 
at  any  lag).  Many  other  models  were  fitted  to  the  data 
with  varying  degrees  of  success  which  were  all  less  than 
that  achieved  by  the  two  models  listed.  The  column 
entitled  "S.E."  represents  the  forecast  standard  error 
for  the  given  model.  In  Table  I,  two  such  values  are 
shown  per  model.  The  reason  for  this  is  to  illustrate 
that  a  given  model  can  achieve  varying  amounts  of 
effectiveness  at  predicting  the  series  values  depending 
upon  where  the  initial  point  of  forecast  takes  place. 

In  the  case  of  the  100  day  model,  two  different  requests 
for  forecasting  (predictions)  were  made.  The  one  which 
resulted  in  a  S.E.  of  157.49  began  at  the  80th  day  of 
the  100  days  of  data  while  the  other  which  yielded  an 
S.E.  of  174.29  began  at  the  71st  day.  The  forecasts 
were  requested  beginning  at  these  particular  days  in 
order  to  allow  the  available  SEESSD  series  values  to  be 
compared  to  the  values  predicted  by  the  models. 
Similarly,  the  365  day  model  resulted  in  an  S.E.  of 
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227.34  and  241.04  for  requested  forecasts  beginning  at 


i 


,'Vv 

.v*V 


the  350t.h  and  200th  days  of  the  data  respectively.  The 
approximate  serif's  mean  for  each  of  the  channels  is 
shown  at  the  bottom  of  each  of  Tables  I  -  IV.  This  is 
shown  in  order  to  give  the  reader  an  idea  of  the 
magnitude  of  the  forecast  standard  error  as  compared  to 
the  mean  series  value.  These  means  are  only  approximate 
since  they  were  initially  determined  over  the  entire  4 
years  of  data  originally  presented  for  analysis  by  LANE. 
Thus,  the  mean  of  each  series  will  change  depending  on 
the  amount  of  data  (number  of  days)  to  which  a  model  is 
fit.  It  was  assumed  here  that  the  means  of  the  lesser 
series  do  not  vary  substantially  from  those  determined 
for  the  entire  four  years  of  the  data.  This  is  a 
reasonable  assumption  since  each  of  the  series  for 
SEESSD,  SEEI,  SEEII,  and  SEEIII  appear  to  be  stationary 
(ie.,  the  mean  does  not  vary  considerably  over  time). 

The  models  shown  in  Tables  II  -  IV  are  displayed 
in  a  similar  manner  to  those  in  Table  I  and  thus 
represent  the  best  models  for  the  SEEI,  SEEII,  and 
SEEIII  channels.  As  i n  the  case  of  the  SFESSD  channel, 
each  of  the  models  for  the  other  three  channels  appears 
to  be  best  represented  by  some  type  of  ARMA  model. 
Indeed,  a  closer  look  reveals  that  almost  without, 
exception  each  of  the  models  for  the  different  channels 
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contains  autoregressive  terms  of  first  and  27th  ordci 


along  with  a  moving  average  term  of  first  order.  In 
each  of  the  tables,  the  term  Zt  should  be  taken  to 


TABLE  II 


Suggested  Models  for  the  SEEI  Channel  (3.4  -  4.9  MeV) 


Data  Length 


100  days 
( Jan  86-May  86 ) 


365  days 

(May  8 5 -May  86) 


Suggested  Model 

(1  -  . 7  9  7  5B  -  . 09  7  3B  2  6 

-  .  04 2 1 B  2  7 ) z  t 

=  (1  +  . 0959B ) a  t 

(1  -  . 6  24  3B  -  . 0806B  2  6 

-  . 0844B  2  7 ) z  t 

=  (  1  +  . 4  6  2  3B ) a  t 


S  .E  . 
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1  .  598 


Approximate  series  mean  =  1.05 


TABLE  III 


Suggested  Models  for  the  SEEII  Channel  (4.9  -  6.6  MeV) 


Data  Length 


Suggested  Model 


TABLE  IV 


Suggested  Models 

for  the  SEE  111 

Channe 1  (6.6 

-  9.7  MeV) 

Data  Lenr 

;th 

Suggested  Model 

S  .  E  . 

100  days 
(Jan  8 6 -May 

86  ) 

(1  -  . 9932B  - 
=  (1 

. 1 222B2  7 ) z  t 
+  . 1 307B ) a  t 

.0167 

365  days 
(May  8 5 -May 

86  ) 

( 1  -  . 9850B  - 
=  (1 

. 1 701B2  7 ) z  t 
+  .  1  227B ) a  t 

.0204 

Approximate  series  mean  =  14.01 


mean  the  current  value  of  the  stated  energetic  electron 
channel  as  it  deviates  from  the  overall  series  mean. 
This  term  was  explained  in  the  previous  chapter. 

Trans f e r  Function  Models 

The  '-esult.s  of  the  transfer  function  modeling 
are  presented  in  Tables  V  -  XVII.  Tables  V  -  VIII 
contain  models  suggested  for  each  of  the  energetic 
electron  channels  using  the  entire  extent  of  the 
available  cross  correlated  data  (706  days  or  from  8  May 
83  to  12  Apr  85).  The  models  include  those  for  which 
the  solar  wind  as  well  as  the  IMF  Bz  component  is 
considered  as  the  input  series.  The  output  series  is 
always  considered  to  be  the  stated  channel.  Thus,  in 
each  of  the  suggested  models  in  ti  ^-se  tables,  the  term 
y t  denotes  the  outpu*  series  current  value  (the  channel 
and  the  term  x  t  denotes  the  Input  series  (either  the  B  ? 


component  or'  the  solar  wind).  The  term  associated  with 
at  represents  the  noise  component,  in  the  model.  Tables 
IX  -  XII  are  duplicates  of  V  -  VIII  except  that  the 
models  shown  only  represent  data  covering  the  last  200 
days  of  cross  correlation  (25  Sep  84  -  12  Apr-  85). 

Tables  XIII  -  XVI  are  once  again  the  same  except  that 
they  only  represent  some  70  days  worth  of  cross 
correlated  data.  Tables  XVII  and  XVIII  are  special  in 
that  they  identify  transfer  function  models  for  the 
SEEIV  channel  based  on  a  solar  wind  input.  As  stated 
earlier  in  this  chapter,  an  individual  series  model  for 
the  SEEIV  channel  was  not  identifiable.  Each  of  the 
tables  exhibits  the  r,  s,  and  b  values  used  in  the 
original  difference  equation  along  with  the  standard 
error  of  the  forecasts.  The  r,  s,  and  b  values  were 
explained  in  Chapter  4.  As  is  evident  from  the  tables, 
many  of  the  models  contain  delay  parameters  equal  to  one 
day,  though  some  have  delays  which  differ  significantly 
from  one. 

Table  V  may  be  used  in  illustration.  The  models 
for  the  Bz  and  solar  wind  inputs  are  as  shown .  The  r, 
s,  and  b  used  initially  were  (0,0,1)  and  (1,2,1) 
respectively.  The-  forecast  standard  errors  for  the 
models  shown  are  186.91  and  187.36. 


TABLE  V 


Transfer  Function  Models  for  the  SEESSD  Channel 
TOG  days  (8  May  83-12  Apr  85) 


Input  Suggested  Model 

Bz  y  t  =  - 1 4 . 96Bx  t 

+  ( l/(  . 5434B  +  . 1 9 1 1 B  2 1 ) )a, 

(r,s,b)  =  (0,0,1) 

S.E.  =  186.91 

Solar 

Wind  y,  =  ((.3342B  +  .1372B2  -  .2040B3) 

/(I  +  . 998  IB ) }x i 
+  (  1  /(1  -  . 5  5  78B  -  . 1 8  7  8B  2  7 ) ) a 

( r , s , b )  =  (1,2,1) 

S.E.  =  187.36 


TABLE  VI 


Transfer  Function  Models  for  the  SEEI  Channel 
706  days  (8  May  83  -  12  Apr  85) 


Input  Suggested  Model 

Bz  yt  =  - . 083Bx  t  +  (  1/(1  -  . 5589B  -  .1431B3 

-  . 1 09 1 B  2  4  -  . 2387B  2  7 ) ) a  t 

( r , s , b )  =  (0,0,1) 

S.E.  =  1.307 


Solar 

Wind  yt  =  (  (  . 0020B 2  -  .0013B4  )/(  1  +  .9998B))xt 
+  (1/(1  -  . 5680B  -  . 04  4  4B 2  -  .1951B3 
+  . 1 5 1 5B  4  -  .  1  906B  2 7  )  la. 

( r , s , b )  =  (1,3,1) 

S.E.  =  1.497 


72 


TABLE  VII 


Transfer  Function  Models  for  the  SEEII  Channel 
706  days  (8  May  83  -  12  Apr  85) 


Suggested  Model 

0088Bx  t  +  (  1/(  1  -  . 6468B  -  .1165B3 

-  . 1 9 SOB 2 ' ) )a  i 

( r , s , b )  =  (0,0,1) 

S.E.  =  .1699 


Solar 
W i nd  y  t 


0002B  2x  t  +  (  1/(  1  -  .  6974B  +  .1173B2 
-  . 1 039B  3  -  .  1 604B  2  7 )  )a, 


( r , s , b )  =  (0,0,2) 
S.E.  =  .1689 


TABLE  VIII 


Transfer  Function  Models  for  the  SEEIII  Channel 
706  days  (8  May  83  -  12  Apr  85) 


Suggested  Model 


B  z  yt  =  ( - . 0  7 1 3B  5  +  ,1351B6)xt 

+  (  1/(  1  -  . 34  34B  - . 1 77 1 B  2  -  .1551B3 
+  . 1 9 3 5B  4 )  )at 

( r , s , b )  =  (1,1,5) 

S.E.  =  2 . 058 

Solar 

Wind  yt  =  ( ( - . 002  3B 7  -  .0023B*)/(  1  -  .6482B))xt 
+  (1/(1  -  .  3  3  4  3B  -  .  1 8 1 3B  2  -  .1  504B3 

+  .  191  IB4 )  )a t 

(r,s,b)  =  (1,1,7) 

S.E.  =  2. 196 


t'.vi.v'iVy.WMvWA’li.W 
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TABLE  IX 


Transfer  Function  Models  for  the  SEESSD  Channel 
200  days  (2f>  Sep  8-1  -  12  Apr-  85) 


Suggested  Model 


yt  =  -24-OlBx,  +(1/(1  -  . 6 1 25B 
-  .  1973B2  7  Mat 

(r,s,b)  =  (0,0,1) 
S.E.  =  189. 13 


Solar 

Wind 


( (  . 6  7  4  9B  -  . 4985B  3 ) / ( 1  -  .9107B))xt 
+  (  1/(  1  -  . 57  35B ) ) a  t 


(r,s,b)  =  (1,2,1) 
S.E.  =  201.25 


TABLE  X 


Transfer  Function  Models  for  the  SEEI  Channel 
200  days  (25  Sep  84  -  12  Apr  85) 


Suggested  Model 


Bz  yt  =  - . 207  8Bx  t  +  (  1/(1  -  .  4547B  -  .0881B3 

-  . 0876B 2  1  -  .3046B27))at 


( r , s , b )  =  (1,2,1) 

S.E.  =  1.556 


Solar 

Wind  yt  =  ((.0026B  -  .0047B3)/(1  -  .9334B))x, 
+  (1/(1  -  .  4937B  -  .  2  7  4  5B  3  +.2193B* 

-  .  1  864B  2  7 ) a  t 

( r , s , b )  =  (1,3,1) 

S.E.  =  1 .546 


TABLE  XI 


■ 


Transfer  Function  Models  for  the  SF.EII  Channel 
200  days  (25  Sep  84  -  12  Apr  85) 


Input 


Solar 

Wind  y » 


Suggested  Model 

0 1 1 5Bx  t  +  (1/(1  -  . 5 7 1 8B  -  .3042B3 
+  . 2504B  4  -  . 306 5B  2  6 ) ) a  t 

(r,s,b)  =  (0,0,1 ) 

S.E.  =  .1889 


0006B  2  x  t  +(  1/(1  -  .  506  IB  -  .2915B3 
+  . 1 598B  4  -  .  344  5B  2  7 ) ) a  t 


(r,s,b)  =  (0,0,2) 


TABLE  XII 


) 


Transfer  Function  Models  for  the  SEEIII  Channel 
200  days  (25  Sep  84  -  12  Apr  85) 


Input 


Suggested  Model 


B  z  Parameter  estimation  terminated,  no 

apparent  model. 

Solar 

Wind  y,  =  .0022B9Xt  +  (1/(1  -  .8535B  +  .1603B2 
-  .  04  39B  6  +  . 3 1 2  5B  1  3  -  .4  784B14))at 

( r , s , b )  =  (0,0,9) 

S.E.  =  1.221 


ri .  fH  Ji. 


TABLE  XIII 


Transfer  Function  Models  for  the  SEESSD  Channel 
Bz  =  70  days  (30  Oct  83-7  Jan  84) 

Solar  Wind  =  71  days  (1  Feb  85  -  12  Apr  85) 


Input  Suggested  Model 

Bz  No  significant  cross  correlations,  no 

apparent  model. 

Solar 

Wind  yt  =  .2161B3xt  +  (1/(1  -  .799lB))at 

( r , s , b )  =  (0,0,3) 

S.E.  =  163.42 


TABLE  XIV 


Transfer  Function  Models  for  the  SEEI  Channel 
Bz  =  70  days  (30  Oct  83-7  Jan  84) 

Solar  Wind  =  71  days  (1  Feb  85  -  12  Apr  85) 


Input  Suggested  Model 

Bz  No  significant  cross  correlations,  no 

apparent  model . 

Solar 

Wind  yt  =  .0005B2xt  +  (1/(1  -  .7357B 

+  . 2 1 1 4B  2 ) ) a  t 


(  r , s , b )  =  (0,0,2) 

S.E.  =  1 .218 


TABLE  XV 


Transfer  Function  Models  for  the  SEEII  Channel 
Bz  =  70  days  (30  Oct  83-7  Jan  84) 

Solar  Wind  =  71  days  (1  Feb  85  -  12  Apr  85) 


m 


Suggested  Model 


(  - • 0 1 9 1 B  1  1  +  . 04  5  3B 1 2 ) x  t 
+  ( l/( . G368B  - 


1 9 1 3B  4 )  )  a 


Solar 
W i nd  y  t 


( r , s , b )  =  (0,1,11) 

S.E.  =  .1244 


( (  . 0007B  3  -  . 00 1  OB  6  )  / (  1  +  .  337GB) ) x  t 
+  (1/(1  -  . 5547B) ) a , 


( r , s , b )  =  (1,4,2) 


TABLE  XVI 


Transfer  Function  Models  for  the  SEEIII  Channel 
Bz  =  70  days  (30  Oct  83-7  Jan  84) 

Solar  Wind  =  71  days  (1  Feb  85  -  12  Apr  85) 


Input  Suggested  Model 

Bz  No  significant  cross  correlations,  no 

apparent  model. 

Solar 

Wind  y,  =  (.00009B2  +  .00005B3  +  .00003B4)xt 

+  (  1  /(1  -  . 3589B  +  .  27 1 5B 2 
-  .  4446B2  5  -  .  4G00B2 1 ) )a , 

( r , s , b )  =  (0,2,2) 

G 


TABLE  XVII 


Additional  data  concerning  the  autocorrelation 


and  partial  autocorrelation  functions  of  the  two  input 


scries  (the  IMF  Bz  and  the  solar  wind  speed)  are 


contained  in  Appendix  G.  These  ACFs  and  PACFs  were  used 


to  identify  the  input  series  necessary  for  the  initial 


part  of  the  transfer  function  modeling  as  accomplished 


by  the  prewhitening  method.  Both  of  these  input  series 


were  shown  to  be  AR ( 1 )  processes.  Also,  Appendix  H  is  a 


listing  of  the  70  day  cross  correlation  functions 


between  these  input  series  and  the  various  energetic 


electron  channels  (output  series)  so  that  the  reader  may 


better  grasp  the  procedures  for  identifying  the  r,  s, 


and  b  values  in  transfer  function  modeling 


Discuss i on 


Apart  from  the  time  series  models  developed  in 


this  study,  good  support  for  some  of  the  results 


obtained  in  prior  studies  of  energetic  electron  flux  is 


evident.  Perhaps  the  two  best  examples  of  this  are  (1) 


the  fact  that  a  large  positive  autocorrelation  showed  up 


very  consistently  near  lag  k  =  27  in  the  ACFs  of  the 


lowest  three  energetic  electron  channels  studied 


(SEESSD,  SEEI ,  SEEII),  and  (2)  the  delay  factor  for  many 


of  the  transfer  function  models  was  on  the  older  of  one 


to  two  days.  Finding  (1)  above  is  in  good  agreement 
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with  McCormick  who  stated  in  his  conclusions  that 
"whatever  processes  affect  electron  f  1  uc  tua  t.  i  ons  in  the 
three  lowest  channels  are  different  from  those  affecting 
the  two  highest  channels"  (McCormick,  1 98  1  : G2 , C 3  )  .  In 
addition,  it  supports  the  findings  of  Paulikas  and  Blake 
who  noted  the  variability  of  energetic  electron  flexes 
associated  with  the  27  day  solar  rotation  period 
Paulikas  and  Blake,  1978:22).  Finding  (2)  is  also  in 
harmony  with  McCormick’s  finding  that  the  electron 
fluxes  in  the  SEESSD,  SEEI ,  and  SEEII  channels  showed  a 
weak  but  evident  correlation  with  one  and  a  half  to 
two-dav  old  solar  wind  speed  data.  (McCormick,  1984:62). 
Moreover,  Paulikas  and  Blake  are  in  agreement  with  this 
delay  factor  (Paulikas  and  Blake,  1978:11).  Most  of 
these  previous  findings  were  mentioned  in  Chapter  II  of 
this  report. 

Referring  to  the  individual  series  models, 

Tables  I  -  IV  seem  to  suggest  mixed  amounts  of  success 
in  using  time  series  analysis.  With  regard  to  Table  I, 
the  best  models  identified  yielded  a  forecast  standard 
error  which  is  undoubtedly  too  great  to  make  these 
models  of  any  practical  use.  A  95%  confidence  interval 
for  a  prediction  one  day  in  the  future  using  either  of 
the  two  models  listed  would  be  an  extremely  broad 
interval  of  about  3  to  4  times  the  mean  value  of  the 


series  itself.  One  should  bear  in  mind  that  the 
standard  error-  of  the  forecasts  depends  somewhat  upon 
where  (ie. ,  beginning  on  what  day)  the  requested 
forecast  is  made  as  is  illustrated  by  the  two  different 
standard  errors  listed  for  each  model  in  Table  I. 

Still,  this  appears  to  render  the  models  unacceptable. 

A  more  reasonable  model  would  give  a  95%  Cl  of,  say,  no 
more  than  half  of  the  mean  value  of  the  series.  With 
this  in  mind,  it  then  appears  that  the  models  in  Tables 
II  and  III  are  likewise  of  little  value. 

The  results  shown  in  Table  IV,  however,  offer 
some  hope.  As  shown,  the  models  for  the  SEEIII  channel 
both  yielded  very  small  forecast  standard  errors  (.0167 
and  .0204  for  the  100  day  and  the  3C5  day  models 
respectively)  as  compared  to  their  respective  series 
means.  This  would  appear  to  make  them  of  possible 
future  use,  but  one  must  keep  in  mind  the  data  as  it 
appeared  in  this  channel.  The  SEEIII  values  were 
characterized  by  fairly  constant  count  rates  (usually 
about  0.12  -  0.15)  except  on  a  few  rare  occasions  where 
the  count  rate  skyrocketed  to  values  such  as  1417.89 
which  was  the  highest  value  recorded  in  this  series. 
Outliers  such  as  this  can  have  a  marked  effect  on  the 
series  mean  and  standard  deviation  and  thus  make  a  given 
model  seem  more  appropriate  than  it  may  in  fact  be.  In 


mm 
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truth,  with  the  exception  of  the  few  outliers  in  the 
SFFTIT  channel,  one  might  very  well  be  able  to  forecast 
the  next  two  days  of  count  rates  for  SEE! II  simply  by 
"eyeballing"  the  past  data  and  noting  that,  for  the  most 
part,  it  varies  little  from  day  to  day. 

Many  of  the  transfer  function  models  are 
likewise  of  doubtful  utility.  In  fact,  it  seems  that 
the  standard  error  for  the  forecasts  based  on  any  given 
transfer  function  model  was  often  times  as  much  or 
greater  than  that  of  the  individual  series  models. 

Again,  bearing  in  mind  the  series  means  and  the  forecast 
standard  errors  indicated,  it  appears  that  the  706  day 
models  shown  in  Tables  V,  VI,  and  VII  are  not  useful. 
Once  again  though,  the  models  for  the  SEEIII  channel 
(Table  VIII)  have  acceptable  standard  errors  as  compared 
to  the  series  mean  and  thus  may  be  of  some  possible 
value  in  prediction. 

The  exact  same  set  of  statements  may  be  said  in 
regard  to  the  200  day  models  (Tables  IX  -  XII).  Once 
again  the  200  day  model  for  the  SEEIII  channel  (based  on 
a  solar  wind  input)  possessed  a  small  standard  error  as 
compared  to  the  mean  value  of  the  SEEIII  channel.  The 
reader  will  note  that  a  transfer  function  model  for  the 
SEEIII  channel  using  the  IMF  B2  as  the  input  was  not 
identifiable.  The  computer  terminated  parameter 


estimation  due  to  estimates  of  the  parameters  "becoming 
too  highly  correlated". 

In  recognition  of  the  fact  that  quite  a  bit,  of 
the  input,  series  data  had  to  be  estimated  (ie.,  the  data 
from  the  NSSDC  on  the  solar  wind  and  IMF  Bz  contained 
many  missing  values),  it  was  felt  that  perhaps  the 
transfer  function  models  were  not  getting  a  fair  chance 
at  success.  Thus,  a  further  examination  of  the  data  was 
done  in  an  effort  to  determine  a  "stretch"  of  erough 
consecutive  days  of  unestimated  data  so  as  to  make  the 
transfer  function  models  less  susceptible  to  any  bias 
caused  as  a  result  of  all  the  estimated  values.  For  the 
IMF  Bz  data,  the  longest  stretch  of  "clean"  data  which 
could  be  cross  correlated  with  the  energetic-  electron 
channels  occurred  during  the  dates  from  30  Oct  83  to  7 
Jan  84  (70  days).  For  the  solar  wind,  the  longest 

stretch  occurred  between  1  Feb  85  to  12  Apr  85  (71 
days).  The  results  of  the  transfer  function  models 
using  these  days’  Bz  and  solar  wind  values  are  given  in 
Tables  XIII  -  XVI.  As  is  evident,  results  using  this 
"unbiased  data"  were  only  marginally  better  than  those 
when  the  data  containing  estimated  values  were  used. 

And,  once  again,  the  transfer  function  model  for  the 
SFFIT1  channel  offered  the  most  possible  usefulness. 

The  reader  may  note  that  in  three  cases,  models  tore  m  t 


identifiable  due  to  lack  of  any  significant  cross 
correlations  between  the  input  series  and  the  stated 
energetic  electron  channel.  In  all  three  cases,  this 
occurred  when  the  B2  component  was  used  as  the  input 
series.  This  again  offers  some  support  for  the  findings 
of  McCormick  who  noted  that  little  correlation  existed 
between  the  IMF  B2  and  the  energetic  electron  fluxes 
(McCormick,  1984:63). 

Tables  XVII  and  XVIII  represent  the-  only 
identified  way  of  modeling  the  SEEIV  channel.  In  both 
cases,  no  model  could  be  identified  using  the  IMF  Bz  as 
the  input.  However,  when  the  solar  wind  was  assumed  as 
the  input,  some  models  of  possible  usefulness  were 
obtained  as  the  standard  error  of  the  forecasts  was 
small  when  compared  to  the  series  mean.  Once  again,  ..he 
same  things  which  were  said  about  the  behavior  of  the 
SEEIII  channel  may  also  be  said  about  the  SEEIV  channel: 
The  series  contains  outliers.  Therefore,  careful  use  of 
these  models  is  in  order. 


VI.  Conclusions  and  Recommendations 


Cone lusions 


I 


This  study  had  some  measurable  "positives"  to 
report.  First,  it  reinforced  previous  findings  that 
changes  in  the  levels  of  energetic  electron  flux  lag 
solar  wind  speed  changes  by  one  to  two  days.  This  was 
brought  out  by  the  fact  that  many  of  the  suggested 
transfer  function  models  had  delay  parameters  (b)  of  1 
or  2  which  indicates  that  the  output  (the  given  electron 
channel)  lags  the  input  (the  solar  wind  or  the  IMF  B - ) 
by  a  one  to  two  day  time  period.  Second,  the  fact  that 
the  autocorrelation  functions  of  the  three  lowest  energy 
channels  (in  addition  to  that  of  the  solar  wind)  showed 
a  strong  positive  autocorrelation  around  the  lag  k  =  27 
definitely  supports  previous  findings  that  the  solar 
rotation  eye  le  pi  ays  a  considerable  part  in  the  makeup 
of  the  solar  wind  plasma  and,  consequently,  in  the 
changes  in  energetic  electron  flux.  A  closer  inspection 
of  the  ACFs  shows  that  a  buildup  of  autocorrelations  at 
lags  just  before  the  k  =  27  day  point  occurs  followed 
generally  by  a  "builddown"  in  the  autocorrelations  after- 
this  point.  A  possible  theory  on  the  reason  for  this  is 
the  surmised  existence  of  a  "hot  spot"  of  the  sun  which 
emits  greater  than  average  amounts  of  plasma  and  which 
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slowly  moves  across  the  solar  surface  while  the  sun 
rotates.  Third,  the  fact  that  the  three  lowest  energs 
channels  displayed  this  2  7  day  lag  trend  while  t  h< 
higher  channels  did  not  gives  support  to  previous 
findings  that  the  lower  energetic  electron  channels 
appear  to  behave  differently  from  the  higher  energj 
channels.  Thus,  the  control  mechanisms  for  the  lower 
energy  channels  seem  to  differ  from  those  of  the  higher 
channe 1 s . 

Of  less  positive  note,  it  seems  that  the  entire 
methodology  of  time  series  analysis  did  not  apply  as 
well  as  had  been  hoped  to  the  data  presented  for  study. 
This  might  lead  one  to  become  disconsolate  with  regard 
to  the  Box  and  Jenkins  method.  However,  some 
qualifications  which  affected  the  analysis  have  to  be 
emphas i zed . 

First,  it  should  be  recognized  that  the  data 
were  anything  but  consistent.  The  LANL  data  on  the 
energetic  electron  count  rates  were  by  far  the  best  with 
regard  to  consistency.  Of  the  data  from  LANL  actually 
used,  only  15  days  worth  of  missing  values  had  to  be 
estimated.  However,  this  was  far  less  than  the  data 
from  the  NSSDC  on  the  solar  wind  and  INF  Bz  which 
contained  stretches  of  days  at  a  time  where  no  values 


were  posted.  The  fact  that  data  were  missing  at  all  can 
only  tend  to  have  adverse  effects  on  the  analysis. 

Second,  there  is  the  iterative  or  "artful" 
nature  of  time  series  analysis  which  must  be  considered. 
The  Box  and  Jenkins  method,  as  presented  in  BMDP, 
requires  that  the  user  interact  with  the  program  in  what 
is  frequently  a  very  long  and  time  consuming  set  of 
processes  before  a  model  is  finally  identified.  The 
user  must  master  the  art  of  examining  the 

autocorrelation  and  partial  autocorrelation  functions  to 
unmask  the  true  behavior  of  a  given  tirrie  series. 
Obviously,  the  more  experience  one  has  at  this  sort  of 
thing,  the  quicker  the  process  of  identifying  the 
appropriate  model  becomes.  The  same  may  be  said  more 
emphatically  of  transfer  function  modeling.  Mastery  of 
the  secrets  held  in  the  cross  correlation  function 
between  a  given  input  and  output  series  is  of  the  utmost 
importance.  Without  it,  an  attempt  to  derive  a  transfer- 
function  model  may  be  doomed  from  the  start.  In  the 
transfer  function  modeling  process  as  presented  in  BMDP, 
the  user  is  required  many  trial  and  error  interactions. 
Initially,  assuming  significant  cross  correlations  exist 
between  the  two  series,  a  guess  at  r,  s,  and  b  may  be 
made  and  the  appropriate  estimates  of  the  difference 
equation  parameters  (all  figured  by  hand)  must  be  input 


into  the  computer.  The  computer  will  further  refine 
these  estimates  or  else  "throw  them  out  entirely".  When 
this  occurs,  the  user  is  forced  to  make  a  new  guess  at 
r,  s,  and  b  and  refigure  new  parameter  estimates.  Also, 
a  great  deal  of  the  time  spent  in  transfer  function 
modeling  consists  of  trying  to  iteratively  determine  the 
noise  portion  of  the  model  (ie.,  the  part  attached  to 
ai).  In  addition,  one  is  reminded  in  Box  and  Jenkins 
that  for  transfer  function  modeling,  there  is  not 
necessarily  one  unique  model  (Box  and  Jenkins, 

1976:387).  Therefore,  despite  the  fact  that  the  model 
one  may  derive  seems  poor  (or  even  good),  the  analysis, 
if  carried  out  a  different  way ,  might  well  yield  more 
favorable  results.  Different  surmisals  of  the  r,  s,  and 
b  parameters  can,  for  instance,  lead  to  different 
models.  In  short,  there  is  often  too  much  left  to  the 
skills  (and  whims)  of  the  user  and  not  enough  left  to 
the  computer,  at  least  as  far  as  the  BMDP  implementation 
is  concerned.  Until  a  more  precise  way  of  performing  a 
time  series  analysis  is  developed  so  as  to  remove  a  lot 
of  the  "art",  this  method  may  well  be  scorned  by  many 
when,  in  fact,  it  may  have  the  potential  of  modeling  a 
given  equispaced  series  of  values  very  accurately. 

Third,  one  must  consider  the  physical  processes 
at  work  in  the  magnetosphere.  As  has  already  been 


brought  out  in  the  literature,  the  int.iract.ions  between 
the  solar  wind,  the  IMF,  and  the  energetic  electron 
fluxes  are  anything  but  trivial.  Therefore  in  this 
respect,  the  results  of  this  study  are  definitely  in 
support  of  this  non- s i mpl i st i c  theory,  since  most  of  the 
transfer  functions  derived  here  seem  to  indicate  that 
the  fluxes  depend  quite  heavily  on  something  other  than 
just  the  solar  wind  and  the  IMF  Bz.  If  the  solar  wind 
and  the  Bz  component  were  the  sole  drivers  of  the 
energetic  electron  fluxes,  then  one  would  suspect  that 
the  forecast  standard  errors  ir  all  the  transfer 
function  models  would,  on  the  whole,  be  considerably 
less  than  what  they  turned  out  to  be.  Put  simply,  there 
is  more  at  work  in  the  magnetosphere  with  regard  to 
energetic  electron  flux  than  is  reflected  here. 

Recommendat i ons 

Notwithstanding  the  somewhat  disappointing 
results  of  the  model  fitting  attempted  in  this  study, 
the  author  continues  to  believe  in  the  Box  and  Jenkins 
method  and  also  that  it  has  utility  with  regard  to  the 
further  study  of  energetic  electron  fluxes  and  the 
magnetosphere  in  general.  If  nothing  else,  even  with 
data  that  contained  many  missing  values,  the  various 
time  series  studied  here  indicate  a  good  capability  to 


be  modeled  by  the  Box  and  Jenkins  method.  This  is 
supported  by  the  fact  that  so  many  of  the 
autocorrelation  and  partial  autocorrelation  functions 
exhibit  what  appear-  to  be  "classic"  time  series  model 
tendencies.  Moreover,  support  is  offered  by  the  models 
themselves.  For  instance,  an  examination  of  BMDP 
Example  #3  in  Appendix  E  shows  that  the  forecasts  mad* 
by  the  model  are  at  least  attempting  to  follow  the 
trends  of  the  actual  values  (page  162).  Thus,  it  is 
recommended  that  time  series  analyses  of  such  data 
continue  to  be  carried  out  in  the  future.  One 
possibility  for  shortening  the  long  iterations  involved 
in  the  BMDP  implementation  of  the  method  is  the 
development  and  use  of  more  sophisticated  software  which 
will  eliminate  some  or  all  of  the  guessing  involved  on 
the  part  of  the  user  and  make  the  method  more  the 
application  of  a  precise  set  of  steps  which  yields  a 
very  accurate  model.  From  consultations  with  advisors 
and  peers,  it  is  the  understanding  that  such  computer 
software  may  be  on  the  horizon.  If  so,  it  would 
definitely  be  to  the  advantage  of  the  person  doing  a 
time  series  analysis  to  use  such  programs. 

With  regard  to  future  data,  it  would  be  most 
helpful  if  it  could  somehow  be  made  more  consistent  and 
required  less  by  way  of  estimation  of  missing  values. 


This  is  truly  an  ideal,  however,  and  the  author  real  ir.'  s 
that  it  is  almost  as  difficult  to  predict  when  a  sat. •  llite 
sensor  is  about  to  malfunction  as  it  is  t  o  prcdi-t 
electron  fluxes.  If  "cleaner"  data  cou  1  d  be  obt  aim'd  ,  t  lie- 
possibility  of  doing  a  time  series  analysis  based  on 
hourly  readings  might  be  an  interesting  undertaking  if  the 
improved  software  spoken  of  above  could  also  be  used.  If 
such  a  computer  program  were  to  yield  an  extremely 
accurate  model  for,  say,  a  transfer  function  between  the 
solar  wind  and  the  SEEIV  channel  utilizing  hourly  data ,  it 
might  be  more  plausible  to  ask  for  (and  rely  on)  a 
prediction  24  -  48  hours  (time  periods)  in  the  future. 

The  reader  will  recall  that  usually  such  predictions  carr\ 
with  them  a  very  high  standard  error. 

Additionally,  nothing  can  help  the  potential  user 
of  the  Box  and  Jenkins  method  more  than  a  solid  foundation 
in  the  material  as  presented  in  the  text.  As  a  minimum,  a 
prior  graduate  level  course  in  time  series  analysis  is 
recommended  . 

Finally,  with  regard  to  the  models  derived  here, 
only  those  with  an  acceptably  small  forecast  standard 
error  should  be  considered  for  use.  The  reader  is 
reminded  that,  time  series  modeling  is  not  a  static  process 
in  that  as  new  data  is  recorded,  the  model  should  bi 
updated  to  reflect  the  series’  latest  trends. 
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Appendix  A:  FORTRAN’  Program  to  Read  the  LAND  Da 


program  readom 

This  program  was  written  so  that  the  dut..  on  t  no 
Omni  Tape  could  bo  compressed  into  a  more 
manageable  file  of  daily  averages .  The.-  really 
important  variables  are  the  variables  for  the 
solar  wind  velocity  (v)  and  thosr 

having  to  do  with  the  IMF  ( bxe , bye , bze , bym , bzm ) . 
Note  that  some  variables  have'  been  declared  as 
integers  and  some  as  real  depending  on  how  the 
NSSDC  specified  their  individual  formats  on  the 
Omni  Tape  (eg.,  12  or  F6 . 2 ,  etc.). 

Note  also  that  the  last  20  real  variables 
beginning  with  the  variable  aptsi  are  the  averages 
of  their  respective  variables  which  we  are 
interested  in  (eg.,  aptsi  is  the  average  of  ptsi, 
av  is  the  average  solar  wind  velocity,  etc.). 

integer  flag,yr,day,hr, brn , idimf, ids w, ptsi , ptsp, 
*kp,c9,r,dst,nn(20) , dread 

real  b , f , thb , phb , bxe , bye , bze , bym ,bzm,sigb,sigf, 
*sigbx,sigby,sigbz,t,n,v,phv,thv,sigt,sign,sigv, 
*sigphv,sigthv, 

*  ss ( 20  )  , 

*aptsi,aptsp,akp,ac9,ar,adst,ab,af,athb,aphb, 

*abxe , abye , abze , abym , abzm , at , an , av , aphv , athv 

The  character  variable  junk  was  introduced  as  a 
variable  into  which  any  miscellaneous  data  records 
which  do  not  contain  numerical  characters  could  be 
read  and  discarded.  For  the  data  set  given,  this 
occured  only  once:  on  the  10th  hour  of  the  338th 
day  of  1982,  a  row  of  asterisks  was  present. 

Thus,  this  data  were  discarded  by  having  it  read 
into  the  variable  junk. 

character  junk 

The  general  purpose  directory  tmp  contained  the 
raw  Omni  Tape  data  (hourly  readings  amounting  to  a 
5  megabyte  file).  The  file  name  containing  the 
raw  data  was  file3.  Consequently,  we  must 
open  this  file  in  order  to  manipulate  the  data. 
Also,  the  file  swdat  must  be  opened  in  order  to 
r<=>ad  in  the  newly  averaged  data. 


93 


Thi'  file  swdat  will  be  the  run  pi  rs  s1  d  ami  mo  it 
manageabl  o  file. 

open  (  2  ,  f  i  1  e  -  ’  /  tmp  /  f  i  1  e  3  ’  ) 
upcn  (  3  ,  f  i  1 e  =  ’  swda t  ’  ) 

The  following  are  small  routines  to  sot  the 
summing  variables  back 

to  zero  at  the  beginning  of  a  new  day. 

do  6  kk=l  , 20 

ss ( kk  )  =  0.0 
continue 

do  8  kk-1  ,  20 
nn ( kk  )  =  1 

c  o  ri  t  i  n  u  e 


The  variable  dread  is  a  conditional  termination 
variable  for  the  program. 

dread  =  0 

i f ( dread . eq . 1 )  goto  200 

The  following  are  the  format  statement  and  lead 
statement  that  are  used  to  read  the  raw  dat.a  from 
the  file  f i le3  . 


for  mat.  (  il  ,i2,i3,i2,i4,2i2,2i3,14f6.2,f8.0, 
*f5.1,3f6.1,f8.0,f5.1,3f6.1,i2,il,i4,i5) 

Note  that  the  original  raw  data  on  the  Omni  Tape 
contained  37  different  variables,  many  of  which 
are  superfluous  to  this  study. 

read (2,11) 

♦flag, yr , day , hr , brn , id imf , i dsw , 

*  pts  i  ,pt.sp,b,  f,thb,phb,bxe,  bye  ,  bze  ,  bym  , 

*bzm , s i gb , s i g  f , s i gbx , 

*  s i gby , s i ebz , t , n , \  , ph v , thv , si gt , si gn , s i gv , 

*  s  i  g ph  v  ,  s  i  g  t. h  v  ,  k  p  ,  e 9  ,  r  ,  d  s  t. 

The  averaging  subroutines  aie  called  f».r  tin 
various  variables  which  have  now  been  read  in. 
These  subroutines  simply  add  the 

new  1 y  read  variable  to  the  previous  s  urn  of  the 


sanu  variable  and  calculate  a  nee  a\'ci  age.  Tie 
subroutine  is  called  for  each 

hourly  reading  until  the  last  leading  of  the  da;, 
is  read  in.  At  that  time,  the  dails  axerng-  \ale 
is  comput  ed  and  rt.-ad  into  the  file  sedat  as  t  h< 
average  value  for  that  variable  for  that  da>. 

cal  1  a vg  i  (  pts i  ,nn(  1  )  ,  ss  (  1  )  ,  apt.s  i  ) 
call  avgi  ( pt  sp , nn ( 2 )  , ss ( 2  )  , apt  sp ) 
rail  avgi ( kp, nn ( 3 ) , ss ( 3 ) , akp ) 
on  1 1  avgi ( c9 , nn ( 4 ) , ss ( 4 )  , ac9 ) 
ea 1 1  avgi ( r , nn ( 5 ) , ss ( 5 ) , ar ) 
call  avgi ( dst  , nn ( 6 )  ,  ss ( 6  )  , adst ) 
call  avgr ( b , nn ( 7 ) , ss ( 7 ) , ab ) 
call  avgr ( f , nn ( 8 )  , ss ( 8  )  , af ) 
call  avgr ( thb , nn ( 9 ) ,ss(9) , athb ) 
cal  1  avgr ( phb , nn ( 10) ,ss( 10) ,aphb) 
cal  1  avg  r ( bxe , nn ( 1 1  )  , ss ( 1 1  )  , abxe ) 
call  avgr ( bye ,nn( 12) ,ss( 12) , abye ) 
call  avgr (bze,nn( 13) ,ss( 13) ,abze) 
cal  1  avgr ( bym , nn ( 1 4 )  , ss ( 1 4  )  , abym ) 
call  avgr(bzm,nn(15)  ,ss(15)  ,abzm) 
call  avgr ( t , nn (16) , ss ( 1 6 ) ,at) 
call  avgr ( n , nn ( 1 7 ) , ss ( 1 7 ) , an ) 
call  avgr ( v , nn ( 1 8 ) , ss ( 1 8 ) , av ) 
call  avgr ( phv , nn ( 19) , ss ( 19) , aph v ) 
call  avgr(thv,nn(20) , ss ( 20 ) , athv ) 

The  23rd  hour  of  the  day  is  the  last  reading  o., 
any  given  day.  Thus  at  this  time,  the  data  must 
be  finally  averaged  and  then  written  into  the  fil 
swdat.  A  total  of  22  variables  will  be  written 
into  the  new  file. 

i f ( hr . eq . 2  3 )  then 

f  or mat (i2,lx,i3,lx,2(f5.1,lx),f4.1,lx, 

*  f 3 . 1 , 1 x , 2 ( f 6 . 1  , lx)  , 

*9(f6.2, lx) ,f8.0, lx,f5. 1 , lx,3(f6. 1 , lx) ) 
write(3,15) 

*  vr , day , apts i , aptsp , akp , ac9 , 

*  ar , ads  t , ab ,af,athb,aphb, abxe , abye , 

*  abze , abym , abzm , 

*  at , an , av , aph v , athv 
end  i  f 

The  following  is  the  subroutine  used  to  close  the 
files  when  all  the  data  from  file3  has  finally 
been  read  and  written.  Note  thal  this  routine 
sets  the  termination  variable  dread  equal  to  the 


program  termination  value;  of  1. 


i  f ( hr . eq . 23  . and . day . oq .101. und . y r . oq . 85 )  t  hen 
c 1 ose ( 2  ) 
close) 3  ) 
dread  =  1 
end  L  f 


The  one  data  record  which  contained  miscellaneous 
"junk"  is  identified  with  the  routine  below. 

This  data  is  to  be  discarded.  This  statement  was 
unique  to  the  data  provided. 

i  f ( hr . eq . 10.and.day.eq.338.and.yr.eq.82)  then 
read (2,16)  junk 
f  ormat ( a  1 80 ) 
end  i  f 

A  series  cf  "if"  statements  to  ascertain  if  the 
last  original  raw  data  record  has  been  read. 
Again,  the  values  used  in  these  3  statements  are 
unique  to  the  data  provided. 

i  f ( hr . eq . 23 . and . day . ne . 1 0 1 . and . y r . nc . 8 5 )  goto  5 
i f ( hr . eq . 23 . and . day . eq . 1 01 . and . yr . ne . 85 )  goto  5 
i f ( hr . eq . 23 . and . day . ne . 1 01 . and . yr . eq . 85 )  goto  5 
goto  10 


These  are  the  actual  subroutines  used  for 
averaging.  Note  that  any  variables  which  have  a 
value  of  zero  are  not  used  in  computing  the  final 
daily  average.  This  is  because  zero  was 
designated  by  the  NSSDC  as  the  value  for  a 
variable  when  its  actual  value  was  "missing"  or 
indeterminate . 


subroutine  avgi ( x , nx , sx , ax ' 
real  sx , ax 
integer  x , nx 

if  (x.ne.0)  then 
s  x  =  x  +  s  x 
ax  =  sx/rix 
nx  =  nx  +  1 

end  i  f 
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Appendix  R : 


FORTRAN'  Program  Combining  LAN’I.  and  NSSDC 
Data 


program  combol 

This  program  was  written  to  read  the  values  in  the 
files  formed  from  the  LAN'L  data  (energetic  elec¬ 
tron  flux)  and  the  NSSDC  data  (solar  wind  speed/ 
IMF  values)  in  order  to  form  one  single  time 
synchronized  data  file  for-  the  purposes  of 
performing  transfer  function  analysis  as  described 
in  Box  and  Jenkins.  The  name  of  the  LAN'L  file 
is  "EEDAT3"  while  the  name  used  for  the  NSSDC  file 
was  "swdat" . 


In  the  following  lists  of  variables,  day,  date, 
abzm  (the  average  daily  Bz  component  expressed 
in  Geocentric  Solar  Magnetospher i c  or  GSM  units), 
av  (the  average  daily  solar  wind  velocity  in 
km/sec),  and  SEESSD,  SEEI ,  SEEII,  SEEIII,  SEEIY 
(the  various  channels  of  energetic  electron  flux) 
are  the  variables  of  interest  despite  the  fact 
that  others  were  copied  into  the  new  time  synch¬ 
ronized  file.  The  other  variables  were  either 
superfluous  or  declared  for  convenience  in  execu¬ 
ting  this  program. 

integer  yr , day , grp , date , pts , count , 
rep, do i , poi , i ( 15  ) ,k,jl,j2 
real  apts i , aptsp , akp , ac9 , ar , ads t , ab , af , 

)  athb , aphb , abxe , abye , abze , abym , abzm , 

t  at , an , av , aphv , athv , 

>  GAMIII, 

*  SEEIII , GAMI I ,GAMIV, SEEI , SEEI I , SEEIY, SEESSD 

character  junkl,junk2 

The  variables  count  and  rep  are  counter  variables 
used  to  control  execution  of  this  program.  The 
variables  doi  and  poi  stand  for  "day  of  interest" 
and  "port  of  interest".  The  127th  day  of  the  file 
containing  the  solar  wind  and  IMF  data  (file 
"swdat")  represents  May  8  which  is  the  first  day 
(matched  to  the  file  containing  the  energetic 
electron  flux  data,  "EEDAT3")  of  sw/IMF  values 
which  should  be  copied.  May  8,  1983  is  the  first 

day  of  data  from  file  "EEDAT3” .  File  "swdat"  has 
data  for  days  prior  to  this  so  all  its  data 
prior  to  this  date  is  unecessary  since  no 


c 

c 


cross  correlations  can  be  performed  on  it 
without  data  for  similar-  days  f ron.  fill-  "EEDAT.'i". 
The*  parameter  poi  is  a  variable  to  name  the  port 
through  which  a  file  is  copied  (a  FORTRAN 
particular-  when  writing  a  file  of  data  to  a  ncv. 
file).  Variables  jl  and  j2  are  also  counters. 

count  =  0 
rep  =  1 

doi  =  127 
poi  =  4 
jl  =  0 
j  2  =  0 

Variable  i  is  an  array  used  to  name  the  particu¬ 
lar  dates  in  file  "EEDAT3"  where  missing  values  of 
variables  of  no  interest  to  this  study  will  have 
to  be  read.  Thus,  when  the  program  determines 
that  one  of  these  dates  is  the  date  of  the  data 
being  read,  it  will  read  the  data  according  to 
a  different  format.  In  so  doing,  it  will  use 
the  character  variables  junkl  and  junk2  to  read 
in  the  missing  values  of  these  useless  variables. 
Without  this  caveat  in  the  program,  the  creation 
of  the  new  single  time  synchronized  file  could  not 
proceed.  Notice  that  there  are  15  dates  on  which 
missing  values  of  superfluous  variables  occur. 

The  other  17  dates  where  missing  values  occurred 
were  in  August  1982.  The  reader  will  recall  from 
Chapter  III  that  the  data  up  until  May  b,  1983 
from  LANL  was  discarded  partly  for  this  reason  (so 
many  missing  values  in  one  month).  Also,  discard¬ 
ing  data  up  to  May  8,  1983  left  exactly  three 

years  (1095  days)  which  could  be  analyzed.  This 
was  felt  to  be  more  than  enough  data  for  analysis. 


i  (  1  ) 

= 

830523 

i  (  2  ) 

r 

830524 

i  (  3  ) 

r 

830526 

i  (4  ) 

r 

840216 

i  (5  ) 

r 

840217 

i  (6) 

r 

840218 

i  (  7  ) 

r 

840620 

i  (8) 

= 

840718 

i  (9) 

= 

841015 

i  (  10) 

=  841016 

i  (  1  1  ) 

=  841030 

i  (  12) 

=  841031 

i  (  13  ) 

=  841101 

i  (1  4  ) 
i  (  1  5  ) 


841102 

841103 


The  following  statements  open  the  relevant  files 
File  "CCFLO"  is  the  name  given  to  the  new  single 
file  containing  combined,  time  synchronized 
sw/IMF  and  energetic  electron  flux  data.  The 
reader  will  notice  that  later  on  in  the  program 
two  other  files  ("CCFL1"  and  "CCFL2”)  were  also 
created.  These  were  initially  thought  to  be 
necessary  for  the  analysis  but  in  fact  were  not. 
The  name  "CCFLO"  is  the  author's  acronym  for 
"cross  correlation  file  with  0  days  l_aS"- 
In  other  words,  "CCFLO"  represents  a  file  of 
data  from  both  "swdat"  and  "EEDAT3"  which  is 
time  synchronized  ( ie. ,  has  zero  days  lag). 

open  (2, file  =  ’swdat’) 
open  ( 3 , file  =  ’ EEDAT  3 ’  ) 
open  (4, file  =  ’CCFLO’ ) 

date  =  830508 

read  (2,10)  yr , day , apts i , aptsp , akp , 
ac9 , ar , adst , ab , af , athb , 
aphb , abxe , abye , abze , abym , 
abzm , at , an , av , aphv , athv 

File  "swdat”  is  read  first  since  it  has  all 
the  excess  data  at  the  start  which  must  be 
discarded  (ie.,  all  data  up  until  May  8,  1983). 

format  (i2,lx,i3,lx,2(f5.1,lx),f4.1,lx,f3.1,lx, 

2 ( f 6 . 1 , 1 x ) , 9 ( f 6 . 2  ,  lx  )  , 
f8 .0, lx, f5 . 1 , lx , 3( f6. 1 , lx)  ) 


We  continue  to  read  data  from  "swdat"  but 
not  copy  it  until  May  8,  1983  rolls  around. 

This  is  the  first  day  of  data  from  file  "EEDAT3 
When  we  finally  reach  this  date,  data  from  both 
"swdat"  and  "EEDAT3"  begins  to  be  copied  into 
the  file  "CCFLO". 

if  ( yr . eq . 8 3 . and . day . eq . doi )  then 
count  =  1 

end i  f 


if  (count.lt. 1)  goto  5 


Continue  to  read  data  without  copying  uni  i] 
reaching  the  appropriate  date. 
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The  following  do  loop  takes  care  of  reading 
in  the  appropriate  values  on  those  days  when 
missing  values  for  superfluous  variables  in 
file  "EEDAT3 "  occur. 

As  mentioned  earlier,  without  this  the  program 
could  not  read  all  the  data.  Notice  that,  when 
the  data  is  read  this  way,  character  variables 
junkl  and  junk2  are  assigned  values  for  what 
would  normally  be  the  unused  variables  pts  and 
GAMII  (see  statement  13  below). 

do  12  k  =  1  ,  1 5 

if  ( date . eq . i ( k ) )  then 
j  2  =  k 
k  =  15 

read  (3,11)  grp , date , junk  1 , SEEI 1 1  , 

*  junk2 , SEEI , 

*  SEEII ,SEEIV,SEESSD 
format  ( i 2 , 2x , i6 , 5x , a2 , lx , 

*  f7.4,a7,3(f6.4,lx) ,f8.4) 
end  i  f 

continue 


if  ( j  2 . gt .  j  1  )  goto  1 6 
c 

c  Statement  13  is  the  normal  read  statement  used 

c  to  read  data  from  file  "EEDAT3".  Statement  10 

c  above  is  only  used  on  the  15  particular  dates 

c  where  missing  values  occur, 

c 

13  read  (3,15)  grp , date , pts , SEEI I I , 

*  GAMII , GAMII I ,GAMIV, 

*  SEEI , SEEII , SEEIV,SEESSD 

15  format  ( i 2 , 2x , i 6 , 5x , i 2 , 4x , f 6 . 4 , 

*  3x  ,  f  7 . 4 , 3x  ,  f  7 . 4 , 2x  ,  f  8 . 4 , 4x  , 

*  f 6 . 4 , 4x , f 6 . 4 , 4x , f 6 . 4 , 2x , f 8 . 4 ) 

c 

c  Statement  16  begins  the  sequence  where  the  data 

e  values  from  both  "EEDAT3"  and  "swdat”  arc  actual  1 

c  written  into  the  new  time  synchronized  file, 

c 

16  jl  =  j  2 

if  (count.eq.l)  then 

write  (poi,17)  date , abxe , abye , abze , 


* 

* 


abym , abzm , av , SEESSD , 

SEEI , SEEII , SEEIII , SEEIV 
format  (i6,lx,5(f6.2,lx),f6.1,lx, 

*  f8.4,lx, 2(f6.4,lx) ,f7.4, lx, f6. 4) 

endi  f 

The  following  list  of  statements  are  termination 
statements.  The  last  day  for  which  time  synch¬ 
ronized  data  may  be  obtained  is  April  12,  1985. 
This  is  the  last  day  of  data  occurring  in  file 
"swdat'1  and  thus  represents  the  last  day  for 
which  we  wish  to  write  down  values  into  the  new 
file.  April  12,  1985  corresponds  to  day  101  in’ 

the  creation  of  file  "CCFLO".  It  corresponds  to 
day  100  and  day  99  in  the  creation  of  files 
"CCFL1 "  and  "CCFL2 "  respectively.  These  last  two 
files  were  net  used  in  the  analysis.  "EEDAT3" 
contains  data  for  dates  after  April  12,  1985,  but 

since  "swdat"  does  not,  we  have  copied  as  many 
values  as  we  can  for  the  purposes  uf  transfer 
function  modeling. 

if  ( yr . eq . 8 5 . and . day . eq . 1 0 1 . and . rep . eq . 1 )  then 
count  =  0 
jl  =  0 
j  2  =  0 

close  ( 2 ) 
close  ( 3 ) 
close  ( 4 ) 
rep  =  2 

endi  f 

If  we  have  not  read  the  last  day,  then  rep  re¬ 
mains  equal  to  1  and  we  continue  to  read  by  re¬ 
turning  to  statement  5  to  iterate  the  process. 

If  we  have  read  the  last  value,  then  rep  is  set 
equal  to  2  (see  the  "if"  statement  above),  and 
the  new  file  is  closed.  We  then  proceed  on  to 
create  and  fill  data  files  "CCFL1"  and  "CCFL2" 
in  the  same  manner  as  we  did  "CCFLO". 

if  (rep.eq.l)  goto  5 
if  (count. eq.l)  goto  20 

if  (rep.eq.2)  then 

open  (2, file  =  ’swdat’) 
open  (3, file  =  ’EEDAT3’) 
open  (7, file  =  ’CCFL1’) 
do i  =  126 


if  ( yr . eq . 85 . and . day . eq . 100. and . rep . eq . 2 )  then 
count  =  0 

jl  =  0 
j  2  =  0 

close  ( 2  ) 
close  ( 3  ) 
close  ( 7  ) 
rep  =  3 

endi  f 

if  (rep.eq.2)  goto  5 

if  (count. eq.l)  goto  25 

if  (rep.eq.3)  then 

open  (2, file  =  ’swdat’) 
open  (3, file  =  ’ EEDAT3 ’ ) 
open  (8, file  =  ’CCFL2’) 
doi  =  125 
poi  =  8 

endi  f 

if  ( yr . eq . 85 . and . day . eq . 99 . and . rep . eq . 3 )  then 
close  ( 2 ) 
close  ( 3  ) 
close  ( 8  ) 
rep  =  4 

end  i  f 


if  (rep.eq.3)  goto  5 


Appendix  C:  FORTRAN  Program  to  Write  Out  All  Pertinent 
Data 


??.: GF.A“  r ATSHC 


THIS  PROGRAM  WRITE?  OUT  All.  THE  LATA  FOR  EACH  OF  T-E 
['AIL?  AVERAGE  VALUES  EECCEDEF  FOR: 

SOLAR  WINE,  EZ,  SEESSD  SEEI,  SEEII,  SEHIII ,  ANL 
SEE! V .  THE  BEGINNING  DATE  IS  KAY  8  ,  1983.  NET ITS 
THAT  VALUES  FOR  THE  SOLAR  WIND  AN!  BZ  COMPONENT  ENI 
ON  •£  APR  81.  THIS  706  DAYS  REPRESENTS  THE  EXTENT 
OF  THE  CSOSS-COERELATABLE  DATA  BETWEEN  THE  SOLAE  WIND/ 
EZ  CONSIDERED  AS  INPUTS  AND  THE  INDIVIDUAL  ELECTRON 
FLIT  CHANNELS  CONSIDERED  AS  OUTPUTS  FOR  TRANSFER 
FUNCTION  MODELING. 


INTEGER  DATE,  COUNT 

REAL  ABZM,  AV,  SEESSD,  SEEI,  SEEII,  SEEIII,  SEE!’/ 
CHARACTER  JUKK1 ,  JUNES,  JUNES 


COUNT  :  1 

OPEN  ;2,  FILE  :  ’ CCFLO ’ ! 
OPEN  13,  FILE  :  ’ EBDAT3 ’ ) 


I  READ  (2,10!  DATE,  JUNE’,,  ABZM,  AV,  SEESSD,  SEE!,  SEE!:, 

*see::i ,  seeiv 


10  FORMAT  ( 1 6 , AL9 , EE . 2 , 1 X , P6 .  L ,  1 X ,  F8 . 4 , ! X ,  Z ( F? .  4 , ! I . . F? .  4 ,  II .  FC .  < 


IF  (COUNT. ED. 1 )  THEN 
WRITE  (6,11) 


FGEMAT  (3I,’DATE’,6I,’BZ’,4E,’SWVEL’,3I,’SEESSD’,4I,’SSE;\ 


4X ,  ’SEEII’ ,31, ’ SEEIII’ ,31, ’SEEIV’ ,/) 


COUNT  :  COUNT  *  1 


WRITE  (6,12)  DATE,  AEZM,  AV,  SEESSD,  SEEI,  SEEII,  SEEIII,  SEEIV 


12  FORMA'  (2X , I £ , 2X , FG . 2 , 2X , F€ . I , 2X, F8 . 4 , 2X,F£ . 4 , 21 , F6 . 4 , IX , 
t  F7 . 4 , 2X , F6 . 4 ! 

IF  ( COUNT. GE. 45:  THEN 
COUNT  :  I 
END  IF 


t;  f  ^  i "E . f  Ir . s : r  4 1 2  ■  GATr  5 


sea:  ',15:  junei,  cate 

:r  FORMA'  i A 4 . : C 


V7.I.VCV, 


J 


m  ‘  ~  Qr'l1*!  r>r  IT  ^ 

r  E:A'  juke:,  r ate ,  jvke:,  see:::,  juke:,  see:,  see::, 
:  see::,  sesss: 

IC  F.SMA?  : A4 , :S , EX, Al , 42, F£ . 4 , A34 , 2 ! f 6 . 4 , 42 ; , Fc . 4 , 21 , FS . 4 ! 

;r  ; CCl’ST.ES . !  *  THEN 
KITE  (5,22! 

22  FORMAT  (3X. ’DATE’ ,6X, *BZ’ ,4X, ’SWVEL’ ,31, ’SEESSD’ ,4X, ’SEE 
»  42, ’SEEII’ ,31, 'SEEIII' ,32, ’SEEIV” , /) 

end:f 


CO UK?  :  COUNT  4  1 

WR:TE  (6,25!  DATE,  SEESSD,  SEEI,  SEEII,  SEE: I I ,  SEEIV 
25  FORMAT  ( 21 , 1 6 , 1 8X , P8 . 4 , 21 , F6 . 4 , 2X , F6 . 4 , 21 , F7 . 4 , 21 , F6 . 4 ! 

IF  (COUNT. GE. 45!  THEN 
COUNT  :  1 

en::f 


IF  (DATE. NE. 860506 !  GOTO  17 


a 


.vvyvv 


Append 

ix  D: 

Listing 

of  Ail  Per 

tirient.  Data 

r  •  *r 

c; 

SWVEI 

CCTQC; 

SEE' 

S 7 ; ! 

SEE"; 

2.23 

4"'.? 

2"6.5CL; 

[ .  P  65 

3.2641 

2. 0344 

A  A  A  J*  r 

830559 

*  r  a 

i.  .  v  v 

474.2 

126.5275 

r.4118 

C.C623 

2.7315 

0.2527 

53C510 

464.6 

344. 5345 

1.3250 

A  f  r?  p 

2.0581 

r  Mr; 

S3C51: 

0.59 

518.2 

11.841! 

0.0597 

0.029" 

C .  C "  8  7 

0.1887 

8  3  C  L !  2 

.•Vi. 

671.1 

110. 3736 

0.2289 

0.055’’ 

0.078? 

3.1812 

830513 

-1.02 

671.1 

282.7424 

ins:-: 

fl  1*11 

2.2253 

0.1886 

830314 

-1.02 

671.1 

429.1727 

2.0434 

0.2541 

0.0923 

0.1932 

830515 

-1.02 

671.1 

396.3473 

2.7705 

0,2319 

0.3945 

0.1933 

330516 

-1.02 

671.1 

374.1903 

3.3052 

0.3017 

0.0981 

0.1937 

830517 

-1.02 

495.6 

143.9721 

1.0575 

0.1227 

0.0628 

0.18C2 

830518 

2.68 

450.8 

60.4168 

0.5044 

C.0692 

0.0799 

0.1794 

830519 

2.37 

410.4 

80,2253 

0.5334 

0.0722 

0.0916 

0.1897 

830530 

3.63 

389.4 

187.6782 

2.061? 

0.2474 

0.0968 

0.2010 

82P521 

5.82 

5  5  7.5 

3.6556 

0.0399 

0.0377 

C  .0816 

0.1868 

830522 

-3.05 

591.8 

31.3407 

0.1020 

0.0433 

0.0800 

0.1849 

83C523 

-2.26 

585.2 

80,7859 

0.4593 

0.0511 

0.0826 

0,1925 

820524 

-1.23 

614.2 

184.3976 

0.9290 

0.1586 

28.4798 

C .2352 

830525 

5.95 

648.0 

184.4024 

0.9290 

0.1586 

28.4'8: 

0.2354 

820526 

5.95 

648  .0 

55.8152 

0.7328 

0.1268 

£ .  C  9  3  8 

0.1361 

830527 

5.95 

648.0 

184,4120 

0.329  : 

3.1586 

2:  d?r 

f  n  a  c  c 

820529 

5.95 

648.0 

70.4957 

0.206? 

0.0541 

f  r  5.  C  r 

6. 15?  5 

830525 

5.95 

64  8.0 

152. 820? 

C.E294 

O.C792 

7 .  ■ 945 

-  err, 

520522 

0.66 

291.3 

119.9056 

Cd7”' 

a  c.njfi 

J  .  .  .  ZL 

r  r  f  ■'  r 

f  '  *  *  C 

820521 

0.79 

420.6 

20.3235 

0.0521 

a  r  f  r.  n 
OilUv 

0.7661 

2.1393 

830601 

-0.86 

432.6 

22.8509 

0.0924 

C.C469 

G .  0  ?  3 ! 

5.2046 

830562 

0.82 

460. 7 

22.5162 

0.070S 

0.043! 

0.7346 

M254 

830622 

0.31 

502.0 

34.9108 

0.130? 

9.9518 

7.0515 

0.2022 

830604 

0.5! 

438.9 

45.3721 

0.165? 

0.0523 

0.7680 

0.2732 

820606 

0.89 

440.9 

55.9589 

0.1887 

0.0535 

0.0889 

0.2040 

820606 

2.72 

517.3 

55.6194 

0.1645 

0.0518 

0.0365 

0.1395 

ClAf A" 

5  v  o  C  w  i 

2.72 

517.3 

29.4189 

0.0874 

0.0464 

0.0865 

0.1999 

830603 

2,i2 

51?  3 

30.8373 

0.0855 

0.0465 

0.0356 

C.2025 

820609 

2.72 

517.3 

15.6942 

0.0588 

0.0423 

0.C82! 

0.1938 

830610 

2.72 

517.3 

4.6172 

0.0436 

0.0434 

0.0813 

0.1905 

830611 

2.72 

357.0 

15. pcs? 

0.0478 

0.0422 

0.0029 

0.1924 

830612 

2.72 

365.2 

4.9617 

0.0426 

0.0408 

0.0635 

0.1914 

8206 12 

-2.81 

555.7 

8.6462 

0.0436 

0.0391 

0.2804 

0.1 897 

8206  1  4 

0.84 

462.5 

87.8682 

0.0656 

0.0462 

j\  a  c  7  9 

0.1999 

820615 

0  ro 

“  i.  .  V  V 

4  37  .9 

63.2878 

0.0922 

0.0495 

a  A  a  r,  r 

3.2068 

820616 

f  t 

385.0 

95.5252 

C .  1762 

c .  0  r  0  7 

r  ' 7  n 

0.2250 

82C61" 

0 

429.8 

63.6880 

6.1423 

0.0162 

3.2595 

0.2016 

820616 

0 .  C  5 

601.3 

24.3370 

C.0742 

0.64*4 

7.2 557 

0.2545 

82C£ I 5 

-1.73 

592.4 

73.4669 

a  1 r- 

t  .  J  b  V  V 

0.2462 

0.055? 

0.202  1 

820622 

-1.73 

476.1 

190.6525 

5 .  4GSS 

ft  r\  r  .3  f 

•J  •  V  V  *  *t 

C .  2  9  4  - 

C.20:’ 

W2W 


wtwir1 


ww 


•w 3 


H2.2236 
244.526! 
1 ?£ . 392 • 
1  1 9 . 8  3  9  S 
273.1828 
128.2135 
67.3353 
100.8528 
41.0516 
12.4443 
48.8261 
58.2352 
33.7165 
24.056! 
20.3643 
8.3975 
14.3057 
22.7935 
19.5251 
5  1.3236 
57.6838 
43.1700 
98.2269 
125. 5243 
183.335! 
3.7769 
51.720? 
56.9284 
212.0317 
286.5022 
279.4203 
124.2042 
4.2854 
21.2244 
96.5466 
271.1768 
252.7928 
72.4116 
52.2086 
59.7341 
115.6035 
192. 2115 


0.471? 
0.5665 
0.e379 
0.5659 
1.7312 
C .  3403 
0 . 22 7  3 
0.3518 
0.0925 
0.0545 
0.1217 
0.1195 
0.099? 
0.0522 
0.0629 
0.0534 
0.0544 
0.0577 
0.0672 
C.C929 
0.107! 
0.1060 
0.0975 
0.1024 
0.1814 
0.0461 
0.0723 
0.1151 
0.3193 
0.6035 
0.5114 
0.5023 
0.0436 
0.0676 
0.1123 
0.3819 
0.3332 
0.1126 
0.0873 
0.1151 
0.1639 
0.2423 
0.0460 


C  ,0605 
0.05CE 
0.0625 
0.0453 
0.0436 
0.0471 
0.0472 
0.0458 
0.045? 
0.0465 
3.0444 
0.0455 
0.0445 
0.0447 
0.0512 
0.0505 
0.0434 
0.0515 
0.0516 
0.054? 
0.043? 
0.0479 
0.0471 
0.0575 
0.0700 
0.0692 
0.062? 
0.0416 
0.0437 
0.0489 
0.0628 
0 .0601 
0.0504 
0.0476 
0.0482 
0.052" 

f\  f  «  « 

0.0422 


0.1012 

2.2139 

0.0902 

0.2266 

0.0882 

0.202? 

0.09' 3 

C.2P53 

0.0872 

0.2034 

0.0:"5 

0.203" 

0.2906 

0.2046 

0.091? 

2.2062 

0.08": 

0.206? 

c .  0  8 7  g 

2.2101 

0.02-74 

0.2123 

c.oss: 

C  .2275 

0.0900 

0.2120 

0.0314 

0.2103 

0.0916 

0.209? 

0.094" 

0.2170 

0.0963 

0.2173 

0,0326 

0.2120 

0.095" 

C.2376 

0.0382 

2.2204 

0,0996 

0.2222 

0.0324 

0.2115 

0.0916 

0.2141 

0.0916 

0.2172 

0.0994 

0.2205 

6.1026 

9.2249 

0.1012 

0.2271 

0.0947 

0.2181 

0.0828 

0.1954 

0.0883 

C .2055 

0.0910 

0.2123 

0.1011 

0.2233 

0.099? 

0.2246 

0,093? 

0.2153 

0.0912 

0.2150 

0.02-3? 

2.2:6: 

Q  c  c r 

2.2246 

0.122: 

1.2222 

C.2S3-' 

0.2122 

e:  swvel  seess: 


see:  v 


?2D9'.7 

2.33 

375.1 

175.0748 

820918 

2.23 

375.1 

253.5108 

£2  S3 1 9 

2.23 

375  .1 

27.8209 

830920 

-1.3’ 

489.5 

269.5759 

83092! 

-1.37 

430.5 

698.4154 

830922 

-0.66 

382.2 

356.9766 

830923 

-0.75 

323.7 

629.8430 

83C924 

1.41 

380.6 

225.0540 

830928 

2.41 

580. 0 

11.2239 

83C926 

1.54 

663.2 

112.1104 

830927 

0.59 

625.7 

433.6756 

830928 

0.59 

625.7 

583.2333 

830929 

0.59 

625.7 

280,9175 

830930 

0.59 

525.7 

505.033? 

83100! 

0.59 

625.7 

337.1849 

831002 

0.59 

625.7 

11.2763 

831003 

-2.88 

427.2 

82.0711 

831004 

-1.50 

426.0 

17.6032 

831005 

2.64 

458.5 

25.7346 

831006 

2.31 

459.1 

10.7290 

83100? 

1.83 

551.7 

29.5990 

831008 

1.33 

558  .1 

42.8985 

821035 

1.46 

504.7 

61.6628 

83IC10 

1.96 

408.2 

—-a 

cz> 

— a 

83101! 

1.96 

408.2 

34.5335 

831012 

1.96 

408.2 

82.9764 

0  0  i  u  c 

1.96 

408.2 

5.2966 

831014 

1.36 

408.2 

13.8164 

831016 

-0.90 

409.7 

92.8338 

821016 

-0.80 

385.3 

98.9013 

82101? 

-4.22 

424  .0 

20.5290 

831018 

-3.77 

488,4 

35.9681 

831019 

0.02 

464.5 

344.3532 

831020 

0.52 

389.4 

548.9635 

83102! 

;.55 

450.4 

305.3996 

831022 

0.30 

560.6 

98.1038 

931023 

0.30 

510.4 

134.7237 

831024 

0.30 

510.4 

282.6889 

821025 

0.30 

510.4 

444.035! 

821026 

0.30 

510.4 

545.8427 

S21027 

C.20 

510.4 

617.8356 

S  2 1 C  2 ; 

2.70 

334.8 

234.5567 

S  3 ! r  2  9 

-2.24 

£  nr  « 

24 .7027 

82  1  72  ! 

r  f.  • 

487.2 

8  IOC 

0.391! 

0.C620 

0. 100-7 

0.2306 

0.5594 

C.C722 

0 . 1764 

C.2?"2 

C.CSS4 

0.046" 

0.0522 

0.222  9 

0.686" 

0.0808 

0.1068 

0 . 2364 

1.9262 

0.1425 

0.1265 

0.2565 

1.1473 

0.096  1 

0.1105 

3.2420 

3.551! 

0.2139 

0.1312 

0.2607- 

0.9161 

0.0948 

C.1C9 

0.232? 

0.0600 

0.3469 

0.0526 

0.2209 

0.211? 

0.0551 

0.0981 

0.227? 

0.8814 

0.0887 

0.1153 

0.2502 

1.5415 

0.1149 

0.1240 

0.2582 

0.7566 

0.C80E 

0.1115 

0.2455 

3.1832 

0.2404 

0.1215 

0.260? 

1.3829 

0.1213 

0.1126 

0.2446 

0.0602 

0.0458 

0.0553 

0.2222 

0.1806 

0.0523 

0.0998 

0.2293 

0.076! 

0.3492 

0.0548 

0.2287- 

0.0723 

0.C487 

0.0565 

0.2273 

0.0504 

C.C453 

0.0956 

0.227  1 

0.0572 

0.3489 

0.0359 

0.23  7  1 

0.0664 

0.0451 

0,0551 

0.2238 

5.0803 

0.052" 

0.0582 

0.2293 

0.1250 

0.051? 

0.0:53 

7 .2294 

0.1106 

0.0525 

0.0555 

0.225" 

6.1255 

0.6531 

0.02-33 

r\  0  7  n 

V  '  U  0  b  U 

0.0489 

0.0461 

0.0523 

0.224" 

0.0524 

C.0473 

0.0351 

0.2263 

0.1428 

0.0508 

0.1014 

0.2266 

0.2360 

0.6555 

3.0536 

0.2347 

0.0774 

0.0475 

0.0544 

0.2256 

0.1054 

0,0475 

0.0940 

0.2274 

0.7599 

0.0990 

0.1153 

0.2486 

2.5230 

0.2271 

0.1278 

0.2603 

1.7729 

0.1640 

0.1121 

0.2426 

C.  1764 

0.0533 

0.0556 

0.2328 

0.2478 

0.0580 

0.09"6 

0.2303 

0.6840 

C.0770 

0.1103 

0.2467 

1.0826 

0.1035 

C.  1225 

0.2657 

CJ 

CT) 

o> 

CJ 

0.1544 

0.1270 

0.2667 

2 . 33f  2 

0.2180 

0.1 205 

0.2624 

(\  7  *  »  n 

V  .  '  1  ■ 

0.108' 

o.:::: 

r .2422 

0.04 "6 

r\  “.  r  i 
u .  ; *. 

8.2224 

rt  rtr  •  1  r  '  •  r  »  r-  -  r  r  r  r 


112 


c.r 

411.3 

199.6826 

0.5820 

0.r32 

0.1 101 

6.2468 

1  «r: 

270.3 

i  *>  p  Q  f  r 

0.3452 

c.cei: 

0.1054 

0.243r 

2.62 

6 1 : .  9 

IV  *1  1  /  < 

0.058? 

0.0481 

0.088" 

"  .23r5 

1.29 

615.9 

r  'i  "  7  q  7 

0.0732 

0.0626 

p  -r-ir 

'  .2'59 

1.45 

532.4 

77.59?’ 

0.1185 

0.0534 

0.1043 

0.2425 

1.2! 

532.4 

116.8316 

0.1934 

C.3E;? 

0.111? 

'■  .2432 

2.25 

532.4 

158.952? 

0.2866 

0.0639 

0.1058 

0.2435 

1.23 

532.4 

50.6835 

C.0994 

0.C522 

0.0965 

0.228! 

2.92 

522.4 

23.9970 

0.0580 

0.0468 

0.0919 

0.2238 

1.62 

589.1 

82.2664 

0.1600 

0.051! 

0.0346 

2.2228 

1.10 

574.8 

270.4604 

0.5886 

0.0749 

0.1066 

0.2364 

4.25 

590.3 

107.2401 

0.4689 

0.0648 

C.C936 

3.215C 

3.83 

490.3 

54.8620 

0.1577 

0.0512 

0.0912 

0.2145 

2.74 

415.5 

36.4082 

0.1052 

0.0486 

0.0346 

5.2211 

1.10 

478.8 

110.6407 

0.354? 

0.0607 

0.103! 

0.2293 

1.58 

535.9 

98.6521 

C.  4022 

0.0624 

0.1022 

0.2345 

0.14 

622.8 

166.9456 

0.3532 

0.068  1 

0.1023 

0.2357 

0.53 

735.0 

307.0937 

1.0932 

0.1060 

0.1096 

0.2411 

0.34 

735.0 
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Appendix  E: 


BMDP  Example  Problems 


BMDP  Example  Problem  #1,  Typical  Preliminary 

Series  Analysis  . 

BMDP  Example  Problem  #?,  Possible  Good  Model 

for  a  Series  Under  Study  . 

BMDP  Example  Problem  #3,  Typical  Transfer 

Function  Model  for  Two  Interacting  Series 
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22 

n  f;  r  I 

A  '  A 

<*«.' 

0.021 

*  IX  • 

2' 

n  f  f  i 

J  .  1  *  *5 

*  IX  * 

Zf 

0.082 

4  IXX* 

34 

-0.026 

♦ 

II 

• 

35 

-0.018 

+ 

r 

* 

36 

-0.034 

+ 

X! 

4 

37 

0.033 

i 

II 

* 

38 

-0.035 

* 

II 

4 

39 

-0.032 

A 

XI 

♦ 

40 

0.004 

* 

I 

♦ 

41 

-0.030 

+ 

II 

4- 

42 

-0.042 

* 

II 

• 

41 

0.028 

4 

I 

♦ 

44 

C.0C1 

♦ 

r 

i 

45 

-0.019 

+ 

I 

4 

46 

-0.024 

+ 

T 

i 

4- 

47 

-0.030 

4 

II 

♦ 

48 

-C.054 

4 

II 

* 

49 

0.02! 

♦ 

IX 

4 

50 

-0.02" 

4 

XI 

t 

r;  * 

3.012 

♦ 

I 

♦ 

52 

0.022 

+ 

II 

4 

52 

-0.013 

+ 

I 

4 

54 

0.045 

4 

IX 

4 

55 

-0.009 

♦ 

T 

4 

56 

-0.034 

+ 

II 

4 

57 

-0.008 

4 

T 

4 

58 

0.010 

+ 

I 

* 

59 

-0.039 

+ 

II 

4 

60 

-0.025 

* 

II 

4 

ARIKA  VARIABLE  IS  SEEI. 

AROR  :  ’(1 1,(27}'./ 

THE  COMPONENT  HAS  BEER  ADDED  TC  THE  MODEL 

THE  CURRENT  MODEL  HAS 
OUTPUT  VAF.IAELE  :  SEE! 

INPUT  VARIABLE  :  NOISE 

ESTIMATION  RESIDUAL  :  RES  EE  I .  ‘ 


W1 


wmiiwi  wi^mw 


& 


rfwiwpnpnorw 


-\-;- 


FIRST  CASE  S'l'MEFF.  TO  RE  USER 
LAST  CASE  NUMBER  TC  EL  USER 


E  S  T I  y  AT  ’  S  E'  C''‘KI  IT "  E;  ’  '  rAST  S "  A  -  - ' 

RELATIVE  CHARGE  IN  EESILL't;  ST’?*  3F  Si  CASES  LESS 


TEAR  C.JRCIe 


SUMMARY  OF  THE  MODEL 

OUTPUT  VARIABLE  -  SEE! 
INPUT  VARIABLES  -  NOISE 

VARIABLE  VAR.  TYPE  MEAN 

SEEI  RANDOM 

PARAMETER  VARIABLE  TYPE 

1  SEEI  AR 

2  SEEI  AE 

RESIDUAL  SUM  OF  SQUARES  = 
DEGREES  OF  FREEDOM  : 

RESIDUAL  MEAN  SQUARE  : 


TIME  DIFFERENCES 
1-  265 

,CTCR  ORDER  ESTIMATE  ST.  ERE. 

1  1  0.7656  O.C33S 

2  27  D.I04D  0.0545 

836.986936 

335 

2.496466 


T-RATIO 

22.15 

1.91 


ESTIMATION  EY  BACRCASTING  METHOD 

RELATIVE  CHANGE  IN  RESIDUAL  SUM  OF  SQUARES  LESS  THAN  O.IC2C;-D4 

SUMMARY  OF  THE  MODEL 

OUTPUT  VARIABLE  -  SEEI 
INPUT  VARIABLES  -  NOISE 


VAEIAE1E  IS  EESEEI . 
“A ILAC  IS  60./ 


firs:  case  number  to  be  used  : 

LAST  CASE  NUMBER  TO  BE  USED  : 
NC.  OF  OBS.  AFTER  DIFFERENCING  : 
MEAN  OF  THE  (DIFFERENCED)  SERIES  : 
STANDARD  ERROR  OF  THE  MEAN  : 
T-VALUE  OF  MEAN  (AGAINST  ZERO)  : 


AUTOCORRELATIONS 


1-  12 

.24 

-.15 

-.12 

-.17 

.  i n 

-.03 

-.02 

.24 

.03 

.05  - 

02 

.02 

ST.E. 

.05 

.06 

.02 

.06 

.06 

.06 

.  06 

.06 

.06 

.06 

06 

f  r 

13-  24 

.02 

-.05 

.03 

.02 

0.0 

0.0 

-.06 

a  i 

0.2 

_  r  *  . 

r  •> 
u 

r.  a 

ST.E. 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.36 

26 

.26 

25-  36 

.06 

.06 

-.02 

.01 

-,c: 

r.  ? 

.  U  v 

A  r 

AO 

“  .  V*. 

.02 

-.03  - 

M 

i 

A  ft 

ST.E. 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

A  f 
\i  . 

.26 

37-  48 

-.03 

.04 

-.0! 

0.0 

.03 

-.07 

a  r 
*  •  «  L 

-  .0! 

A  • 

.  t  : 

:: 

ST.E. 

.36 

.06 

.06 

.06 

.06 

r  r 

.oc 

.  06 

.06 

a  * 

A  A 

;  l 

.of 

45-  60 

-.06 

-.02 

-.02 

0.0 

.02 

.02 

A  ft 

pr 

-.03 

-.01  - 

Ar 

ft  ft 
«  l  L 

ST.E. 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

.06 

06 

.06 

3 LOT  OF 

A'JTCCC'EEELATT 

ON; 

-1.0 

0.8 

0.6 

0.4 

•0.2 

0.0 

0.2 

0.4 

0.6 

0.8  1 

A 

I 

0.238 

♦  111*111 

2 

-0.155 

I+XII  ♦ 

3 

-0.124 

mi  * 

4 

-0.172 

bill  t 

5 

-0.097 

+XXI  * 

6 

-0.027 

+  V.  - 

1 

-0.025 

t  II  • 

8 

0.038 

*  II  * 

BMD?  Emple  Prolhi  JL 


CcBputer:  CYEER 

Data  set  used:  Last  two  years  tf  EET AT 
Mode!  fit  atteapted:  AS  1,27  MA  ! 

Good  aodel 

IPAGE  1  8MDP2T 


BMDP2T  -  BOX-JENEINS  TIME  SERIES  ANALYSIS 
EMDF  STATISTICAL  SOFTWARE,  IRC. 

1964  WESTWOOD  3LVD.  SUITE  202 
(213!  <75-5700 

PROGRAM  REVISED  APRIL  1982 
MANUAL  REVISED  -  1981 

COPYRIGHT  (Si  1982  REGENTS  OF  UNIVERSITY  OF  CALIFORNIA 

TO  SEE  REMARES  AND  A  SUMMARY  OF  NEW  FEATURES  FOR 
THIS  PROGRAM,  STATE  NEWS.  IN  THE  PRINT  PARAGRAPH. 

THIS  VERSION  OF  BMD?  HAS  BEEN  CONVERTED  FOR  USE  ON 

W  SOCC  AND  CYBER  SERIES  COMPUTERS  BY 

BMDF  PROJECT,  VOGELBACE  COMPUTING  CENTER 

NORTHWESTERN  UNIVERSITY 

2129  SHERIDAN  ROAD 

EVANSTON,  ILLINOIS  5020! 

(312!  <92-368! 

RELEASED  AUGUST  1983  FOE  FTN5  COMPILERS 
EXECUTED  ON  86/10/11. AT  20.09.13. 


PROGRAM  CONTROL  INFORMATION 

/PROELEM  TITLE  IS  ’GEOSYNCH  ENERGETIC  ELECTRON  DATA’ 
/INFU7  VARIABLES  ARE  II. 

FORMAT  IS  FREE. 

RECLEN:97, 

/VARIABLE  NAMES  ARE  PBRICD , DATE , PTS , SEE! II , 

GAMI, GAMUT, GAMIV, SESI, SEEII. 
SEEIV, SEESSC . 

/TRANSFORM  DELETE  :  1  TO  "30. 

/  ENL 

PROELEM  TITLE  IS 

GEOSYNCH  ENERGETIC  ELECTROS  LATA 

NUMBER  OF  VARIABLES  TO  REAL  IN.  .  .  .  . 

NUMBER  OF  VARIABLES  ADO  EL  BY  TRANSFORMATIONS.  . 


•w 


TOTAL  NO  “PER  OF  VASIA?:.-; .  1! 

Nl’HEER  OF  CAFES  TO  REAP  IS . TC  E SI- 

CASE  LABELING  VARIABLES  . 

HISSING  VALUES  CRETE?!'  BEFORE.  OR  AFTER  TRAN'S.  .  NR  I  TREE. 

E LAKES  ARE . KISSING 

INPUT  UNIT  NUMBER .  I 

8EWINZ  INPUT  UNIT  PRICE  TO  REAPING.  .  DATA.  .  .  NC 
NUMBER  OF  KURDS  OF  DYNAMIC  STORAGE .  3999E 


»mt  TRAN  PARAGRAPH  IS  USER  ***** 

1  PAGE  2  BMDP2T  GBOSYNCH  ENERGETIC  ELECTRON  DATA 

VARIABLES  TO  BE  USED 

1  PERIOD  2  DATE  3  PTS  4  SEEIII  5  GAMI 

6  GAM III  1  GAMIV  5  SEEI  9  SEEII  10  SEEIV 

II  SEESSD 
INPUT  FORMAT  IS 
FREE 

MAXIMUM  LENGTH  DATA  RECORD  IS  9!  CHARACTERS. 

NUMBER  OF  CASES  READ .  1095 

CASES  KITH  USE  SET  TC  NEGATIVE  VALUE  ....  730 

REMAINING  NUMBER  OF  CASES  .  365 

I  PAGE  3  EMDP2T  GEOSYNCH  ENERGETIC  ELECTRON  DATA 


ARIMA  VARIABLE  IS  SEEIII. 

AROR  :  ’ ( 1 ) , ; 2? I ’ . 

MAOR  :  ’(!;’■/ 

THE  COMPONENT  HAS  BEEN  ADDED  TC  THE  MODEL 

THE  CURRENT  MODEL  HAS 
OUTPUT  VARIABLE  :  SEEIII 
INPUT  VARIABLE  :  NOISE 

1PAGE  4  BMDF2T  GEOSYNCH  ENERGETIC  ELECTRON  DATA 


ESTIMATION  RESIDUAL  :  RESEEIII ./ 


ESTIMATION  BY  CONDITIONAL  LEAST  SQUARES  HETHGE 

RELATIVE  CHANGE  IN  RESIDUAL  S'JM  OF  SQUARES  LESS  THAN  . IOOOE-C4 

SUMMARY  OF  THE  MODEL 

OUTPUT  VARIABLE  -  SEEIII 
INPUT  VAP.IAELES  -  NOISE 


VAivTAHLE  VAR.  TYPE  HE  AX  TIKE  [IFFEEFX’E: 
SEE!  1 1  RANDOM  1-  HI 


TE 

?.  VARIABLE 

TYPE  FACTOR 

ORDER  ES'IHA'E 

ST .  ERR. 

T-EA.T10 

1 

SEEIII 

MA  ; 

1  -.1223 

.  0552 

-2.21 

2 

SEEIII 

AF.  1 

1  .955: 

.  0095 

102.8! 

3 

SEEIII 

AF.  2 

27  .1594 

.  0545 

3.1! 

RESIDUAL  SUM  OF  SQUARES  : 
DEGREES  OF  FREEDOM  : 
RESIDUAL  HEAR  SQUARE  r 


ESTIMATION  EY  BACECASTING  METHOD 

RELATIVE  CHANGE  IN  RESIDUAL  S UK  CF  SQUARES  LESS  THAN  . lOOC'E-O' 


SUMMARY  OF  THE  MODEL 

OUTPUT  VARIABLE  -  SEEIII 
INPUT  VARIABLES  -  NOISE 

VARIABLE  VAR.  TYPE  MEAN 


TIME  IIFFEEEN.’E: 


RANDOM 

1-  255 

’TER  VARIABLE 

TYPE 

FACTOR  ORDER  ESTIMATE 

ST.  ERR. 

T-RATIC 

1  SEEIII 

MA 

1  1  -.122? 

.0547 

-2.24 

2  SEEIII 

AR 

I  1  .9850 

.  0045 

218.90 

3  SEEIII 

AR 

2  27  .1701 

.0543 

3.13 

RESIDUAL  SUM  OF  SQUARES  :  .130532  (BACRCASTS  EXCLUDED) 

DEGREES  OF  FREEDOM  :  334 

RESIDUAL  MEAN  SQUARE  :  .000351 

1PAGE  5  BMDP2?  GHOSYNCH  ENERGETIC  ELECTRON  DATA 


/ 


>Ar..AsLt  15 


N'.'KEEr. 

OF  OESFEVAT 1 0 

NS 

: 

“LAN  'F 

mi*  j  ■■  n 

•  R  £  i  *■' 

! FFERE 

KCE1 

;  SEN' 

■  r 

K  W  - 

STAN;  Ar. 

:  ERROR 

OF  TH 

l  KEAN 

: 

T- VALUE 

OF  HEAR  { AGA 

INST 

ZERO 

!  - 

AITQCCS 

RELATION’S 

-.03 

-.16  - 

.0? 

-.15  ■ 

-.13  - 

ST.E. 

.03 

.05 

.05 

.05 

.06 

13-  24 

.12 

-.06  - 

.01 

-.0! 

C.O  - 

ST.E. 

.06 

.06 

.06 

.06 

.06 

25-  36 

08 

.11  - 

.02 

0.0 

.04  - 

ST.E. 

.06 

.06 

.06 

.06 

.06 

.or?:: 
! . esc : 


a  : 

V  1 

-.01 

i -  n 

•  V  1 

A  4  A 

.  u ; 

-.02 

r  * 
-  .  I'M 

06 

.06 

.06 

.66 

.06 

.06 

t\  r 

.  t  t 

04 

.01 

.0! 

.03 

-.03 

-.08 

-.07 

06 

.06 

.06 

.06 

.06 

.06 

.06 

01 

-.07 

c— 

o 

.03 

0.0 

-.04 

-.0: 

Of 

.06 

.06 

.06 

.06 

.06 

.06 

PLOT  OP  AUTOCORRELATIONS 


.CSC1 


«  7 

- .  C22 

*  ?:  * 

2? 

.014 

*  I  t 

2S 

.044 

*  12  * 

30 

f"  r 

*  I  + 

3! 

-.O’l 

♦ni  * 

Oi 

.011 

+  in* 

33 

.052 

*  12  * 

34 

i  I  * 

3! 

-.042 

♦  2!  + 

36 

-.084 

+  22:  * 

1  PAGE 

6  BMDF2T 

CEOSYKCfi  ENERGETIC  ELECT: 

FORECAST  CASES  :  24. 

START  :  242./ 


FORECAST  ON  VARIABLE  SEE!  1 1  FROff  TIRE  PERIOD  24! 


PERIOD 

FORECASTS 

ST.  EEE. 

actua: 

54! 

■!2?34 

•  U'vJ 

.  13541 

346 

.1250? 

.02040 

.12521 

34? 

.12266 

.03756 

.13231 

348 

.12352 

.24256 

•  .  j  v «. ! 

349 

.1215: 

.04561 

.1  3501 

Ov0 

.12201 

A  c r,  o  r 

.  V  V  *■'  V 

.  1 3  ’  9 1 

35: 

.12453 

.06702 

.13421 

352 

.12073 

.06063 

.1  350! 

.12124 

.06394 

.1360! 

354 

.11835 

.06699 

.1405! 

355 

.11581 

.06332 

.1409! 

356 

.11300 

.0724? 

.1424! 

35? 

.11136 

.07494 

.1364! 

358 

.11018 

.07727 

.14211 

359 

.10759 

.07946 

.14191 

360 

,10763 

.08153 

.14211 

361 

.10566 

.08349 

.13341 

352 

.10449 

.08535 

.  1423! 

363 

.10270 

.08"  1 1 

.1474! 

364 

.  IC!’’? 

.055'? 

.244:' 

26! 

‘CAM, 

AM  +  * 

<  O /  6 

365 

f  *1 :  r  c 

ft  ^  *  «  • 

n  p  n 

.02-62 

_  r ; ; r 

2  6  s 

f • ;  f  r  - 

r  •  »  "  r 

stamar: 


:et  cos:: 


js'a: 


IrAGE 


?  B“rr*T  G E C S Y K ? H  EKELZETI ?  EL 


en:  ' 

KEY? EE  IF  INTEGER  WC  F.L 5  CF  STCSAGE  UFEI 
CF"  T  L  MF  USE?  SEC 'ST'S 

i  p  i  f  7  g  g  W  r  p  rt  'p 

P“!?2T  -  BOX- JEKEINS  TIME  SERIES  ANALYS 
EXECUTED  OK  86/10/11. AT  DG.C3.c7, 

PROGRAM  CONTROL  INFORMATION 

NO  NOSE  CONTROL  LANGUAGE. 


PROGRAM  TERMINATED 


BKDP2T  -  BCX- JEKKINS  TIKE  SERIES  ANALYSIS 
BMDF  STATISTICAL  SOFTWARE,  INC. 

1964  WESTWOOD  ELVD.  SUITE  202 
(213)  475-5700 

PROGRAM  REVISE:  APSIL  1982 
MANUAL  REVISED  -  198! 

COPYRIGHT  (CJ  lsts2  REGENTS  OF  UNIVERSITY  OF  CALIF! 

TO  SEE  REHARES  AND  A  SUMMARY  OF  NEW  FEATURES  FOE 
THIS  PROGRAM,  STATE  NEWS.  IN  THE  FEINT  PARAGRAPH. 


THIS  VERSION  OF  BMDF  HAS  BEEN  CONVERTED  FOE  USE  ON 

CDC  5001  AND  CYBER  SERIES  COMPUTERS  BY 

BMDF  PROJECT,  VOGELBACE  COMPUTING  CENTER 

NORTHWESTERN  UNIVERSITY 

2129  SHERIDAN  ROAD 

EVANSTON,  ILLINOIS  Oil! I 

[V.l]  492-258: 

RELEASED  AUGUST  1983  FOR  FTN5  COMPILERS 
EXECUTED  ON  38/10/19. AT  10.50.33. 


PROGRAM  CONTROL  INFORMATION 

/PROBLEM  TITLE  IS  ’GECSYNCH  EE  TIMf  S-RIES, 

TRANSFER  FUNC  FROM  SOLAR  WINE ,  IMF  5 
/INPUT  VARIABLES  ARE  12. 

FORMAT  IS  FREE. 

RECL  :  90. 

/VARIABLE  NAMES  ARE  DATE. AEXE. A.BYE.AELE, 

AEYM,  A.BLM,  AV.SEESSE ,  SEE; , 
SEEII.SEEIi;, SEEDY. 


S'UMFER  Of  VABIABLES  TO  READ  IN .  '1 

HUMBER  OF  VARIABLES  A CL  EL  FV  TF.ANSF LR^ATI OKS .  . 

TOTAL  NUMBER  Or  VARIABLES .  1 2 

HUMBER  OF  OASES  TO  REAL  IK . 'LL  ZNI 

CASE  LAE: CISC  VARIABLES  . 

MISSING  VALUES  CHECKED  BEFORE  OR  AFTER  TEAKS.  .  NEITHER 

BLANKS  ARE . MISSING 

INPUT  UNIT  NUMBER .  I 

REWIND  INPUT  UNIT  PRIOR  TC  READING.  .  DATA.  .  .  NO 

N'JHBER  OF  WORDS  OF  DYNAMIC  STORAGE .  <9959 

mil  TRAN  PARAGRAPH  IS  USED  mu 

IPAGE  2  BMDFDT  GEOSYNGH  EE  TIME  SERIES,  TRANSFER  FIT  FROM  SOLA?  KIND,  IK 


VARIABLES  TO  BE  USED 

1  DATE  2  AP5E  3  A3YE  <  A.EZE  I  AEYM 

f  tfiy  ?  av  t  SEESSI  9  SEE!  ID  SEEII 

I!  SEEIII  12  SEED V 
INPUT  FORMAT  IS 
FREE 

MAXIMUM  LENGTH  DATA.  RECORD  IS  $2  CHAEACTESC . 

NUMBER  OF  CASES  READ .  TC' 

GASES  KITH  USE  SET  TC  NEGATIVE  VALUE  ....  90S 

REMAINING  NUMBER  OF  GASES  .  211 

IPAGE  3  BMDP2T  SEGSYKGH  EE  TIME  SERIES,  TRANSFER  F".N2  FE/M  SOLAR 


WILT,  IM 


A.P.IMA  VARIABLE  IS  A.EZE. 

AROF.  :  ’(!!'• 

CENTERED./ 

THE  COMPONENT  HAS  BEEN  ADDED  TC  THE  MODEL 

THE  CURRENT  MODEL  HAS 
OUTPUT  VARIABLE  :  ABZM 
INPUT  VARIABLE  :  NOISE 

IPAGE  I  EMDF2T  GEOSYNCH  EE  TIME  SERIES,  TEANSFER  FUND  FROM  SOLAR  KIND,  IM 


ESTIMATION  RESIDUAL  :  EE. 


ESTIMATION  BY  CCN'ITICKAI  LEA:T  SCARES  METHOD 


RELATIVE  ’HANOI  IN  EACH  ESTIMATE  LESS  THAN 


SwWW 


FILTER  VARIABLE  IS  SEESSI. 

RESIDUAL  =  RY./ 

RESIDUAL  KEAN  SQUARE  :  4EE33.04773I 

VARIABLE  SEESSI-  I?  FILTERED,  RESULTS  ARE  STORED  IK  VARIABLE  RY 
I  PAGE  6  8HEPI?  CEiSEKCH  EE  TIKE  SERIES,  TEAKS?  Er.  F  .S I  FEIV  SI  LAE  VIK 

CCF  VARIABLES  ARE  SX,  RY. 

MAX  LAG  IS  10./ 


EFFECTIVE  N'JKEEE  OF  CASES  :  139 


CCRRELATIOS  OF  EX 

AM  E: 

CO 

CROSS  CORRELATIONS  OF  RX 

(I)  AND  RY 

C+E! 

!-  :e  -.13  .co  -.os 

.It  .11  .01 

.03  -.01  .30  - 

ST.E.  .0?  ,07  ,r 

r  n  {\n  r.r 

o  •»  r " 

.  V  :  «v  .  v-  : 

''I'll  1ATI1-KS  'F  r.V 

(L  AKI  RX 

i  *  -  r  1  -  r  « 

r  r.  ft  <  ft  f 

Q  rr  *  ft "  "  f  “ 

' r  f'n  r  7 

.  a  r,l 

*  - 

TO  THE  MODEL 


CENTERS :./ 

THE  COMPONENT  HAS  EEEK  ADDED 

THE  CURRENT  MODEL  HAS 
OUTPUT  VARIABLE  :  SEESSD 
INPUT  VARIA5LE  =  NOISE 
I  PAGE  9  BKDP2T  GEOSYNCS  EE  TIME  SERIES,  TRANSFER  FUKT  FROM  SOLAF.  NIKI 

INPEF  VARIABLE  IS  AS2M . 

UPOECEES  :  ’(l)’. 

UPVALUES  :  -24 .77275./ 


THE  COMPONENT  HAS  BEEN  ADDED  TO  THE  MODEL 


THE  CURRENT  MODEL  HAS 
OUTPUT  VARIABLE  :  SEESSD 
INPUT  VARIABLE  :  NOISE  AE2M 

1PAGE  10  BMDP2T  GEOSYNCH  EE  TIME  SERIES,  TRANSFER  FUKT  FROM  SOLAR  WIKI 


ESTIMATION  RESIDUALS  :  RYX. 

METHOD  IS  CLS ./ 


ESTIMATION  BY  CONDITIONAL  LEAST  SQUARES  METHOD 
RELATIVE  CHANGE  IN  RESIDUAL  SUM  OF  SQUARES  LESS  THAN 

SUMMARY  OF  THE  MODEL 

OUTPUT  VARIABLE  --  SEESSD 
INPUT  VARIABLES  -  NOISE  AE2M 

VARIABLE  VAR.  TYPE  MEAN  TIME  DIFFERENCES 


PARAMETER  VAR! API :  TYPE  FA  R  'EI'EF.  ESTIMATE  M .  FEE. 

;  apz“  up  :  :  -g . t: :  ?.G24: 

EES  IP  UAL  S’JK  OF  SJ'JASSS  :  iriC2"3.8?’.327 
DEGREES  OF  FREED U".  :  ! SS 

EES : DUAL  KEAN  S4UASE  :  74294. 342"85 

1  FACE  :i  BMDP2?  GEOSTKCfi  EE  TIME  SERIES ,  TRANSFER  FUME  FROM  SOLAR  KIKE,  IMF  El 
ACF  VARIABLE  IS  F. 71./ 


NUMBER  OF  OBSERVATION'S  : 
MEAN  OF  THE  (DIFFERENCED)  SERIES  r 
STANDARD  ERROR  OF  THE  MEAN  : 
T-VALUE  OF  MEAN  (AGAINST  ZERO)  : 


199 

4.9694 

19.318? 

a  r  *»  r 


CORRELATIONS 


I-  12 

.65 

.42 

.20 

.12  .01 

-.OS 

-.10  -.12 

-  f,?  !>'  IS 
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METHOD  IS  CLS . / 


ESTIMATION  BY  CONDITIONAL  LEAST  SQUARES  HETHOI 


RELATIVE  CHANGE  IN  RESIDUAL  SUM  OF  SQUARES  LESS  THAN  .10C0E-Q4 


SUMMARY  OF  THE  MODEL 

OUTPUT  VARIABLE  -  ABZM 
INPUT  VARIABLES  -  NOISE 
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THE  CURRENT  MODEL  HAS 
OUTPUT  VARIABLE  :  AV 
INPUT  VARIABLE  :  NOISE 


ESTIMATION  RESIDUALS  :  EE. 

METHOD  IS  CLS .  / 


ESTIMATION  BY  CONDITIONAL  LEAST  SQUARES  METHOD 


RELATIVE  CHANGE  IN  EACH  ESTIMATE  LESS  THAN  . I0T0E-03 


SUMMARY  OF  THE  MODEL 

OUTPUT  VARIABLE  -  AV 
INPUT  VARIABLES  -  NOISE 

VARIABLE  VAR.  TYPE  MEAN  TIME  DIFFERENCES 
AV  RANDOM  REMOVED  1-  ?I 


PARAMETER  VARIABLE  TYPE  FACTOR  ORCEE  ESTIMATE 
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FILTER  VARIABLE  IS  SEFI. 

RESIDUALS  =  EY . / 

RESIDUAL  MEAN'  SQUARE  : 

VAF.IABiE  SEEI  IS  FILTH F.EI , 


VARIABLES  ARE  EX,  BY. 
MAXLAG  IS  15./ 


EFFECTIVE  NUMBER  OF  CASES  :  70 
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ST.E .  .12  .12  .12  .12  . ' 2  .11  .11  , V  , '  I  .1; 

13-  15  -.06  .04  -.12 

ST.E.  .13  .13  ,13 
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'37.7764’ 
416.53924 
-210  fi'224 
-235.22776 
213,51524 
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The  results  from  this  study  indicate  that  the 
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Finally,  most  of  the  transfer  function  models  derived  with 
the  solar  wind  considered  as  the  input  to  a  given  channel 
resulted  in  delay  parameters  of  about  2  days  between  the 
input  change  in  solar  wind  velocity  and  the  observed  outjiut 
change  in  election  flux  which  support  s>  findings  from  prior 
stud  i  es- .  »  'V-  a 


